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Abstract 

Motivated by recent experiments in ultracold atomic gases that explore the nonequi- 
librium dynamics of interacting quantum many-body systems, we investigate the 
nonequilibrium properties of a Fermi liquid. We apply an interaction quench within 
the Fermi liquid phase of the Hubbard model by switching on a weak interaction 
suddenly; then we follow the real-time dynamics of the momentum distribution by a 
systematic expansion in the interaction strength based on the flow equation method 
[1]. In this paper we derive our main results, namely the applicability of a quasipar- 
ticle description, the observation of a new type of quasi-stationary nonequilibrium 
Fermi liquid like state and a delayed thermalization of the momentum distribu- 
tion. We explain the physical origin of the delayed relaxation as a consequence of 
phase space constraints in fermionic many-body systems. This brings about a close 
relation to similar behavior of one-particle systems which we illustrate by a discus- 
sion of the squeezed oscillator; we generalize to an extended class of systems with 
discrete energy spectra and point out the generic character of the nonequilibrium 
Fermi liquid results for weak interaction quenches. Both for discrete and continu- 
ous systems we observe that particular nonequilibrium expectation values are twice 
as large as their corresponding analogues in equilibrium. For a Fermi liquid, this 
shows up as an increased correlation-induced reduction of the quasiparticle residue 
in nonequilibrium. 
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1 Introduction 



1.1 Experimental approaches 

Recent experimental progress has stimulated the investigation of interacting 
many-particle quantum systems in nonequilibrium. In 1998 ultracold atoms 
confined by a harmonic trapping potential have first been loaded [2] onto 
arrays of standing light waves known as optical lattices [3]. Since that time 
the implementation of paradigmatic quantum models of condensed matter 
physics, foremost the Hubbard model jl] , in such systems with long coherence 
times and high tunability became a feasible task. While in solids parameters 
like the lattice spacing or the interaction strength are fixed constants prede- 
fined by nature, optical lattices provide a technique to experimentally study 
consequences of their time-dependent variation [5J. The observation of the 
equilibrium Mott-Hubbard phase transition in the 3d bosonic Hubbard model 
[6] was soon followed by the discovery of collapse and revival phenomena when 
the system was quenched between its phases, i.e. when sudden changes were 
applied to the particle interaction [7J. Similar nonequilibrium behavior was 
found for quenched Id hard-core bosons jH] which continued to oscillate with- 
out relaxation. 

Experiments with ultracold fermions had to overcome more technical difficul- 
ties. The observation of a Fermi surface in optical lattices [9], of fermionic 
correlations [10], of superfluidity [TT] and the study of interaction-controlled 
transport by a quench in the trapping potential [T2] prepared the recent ob- 
servation of the Mott-Hubbard transition in the repulsive fermionic Hubbard 
model [13] . This rapid progress of experimental sophistication gives hope that 
the predictions presented here will be subject to experimental observation in 
the near future. 



1.2 Thermalization debate 

Such new experimental opportunities and pioneering results have provoked 
further theoretical investigations into the nonequilibrium dynamics of well- 
established many-body model systems, e.g. for spin models |14|I15|I16|I17II18|I19] . 
BCS systems [2THI21|I22||2"3"] . ID hard-core bosons [!M1I25|I26] . a Luttinger liquid 
[271128] . the Richardson model [29], and the Falicov-Kimball model |30j. Ex- 
cited by a quantum quench, i.e. by a sudden switch-on of the interaction term 
in the Hamiltonian, each of them exhibits its individual dynamics; however, 
they share the common feature that in many cases thermalization, i.e. the 
relaxation of time-averaged quantities towards a thermal ground state, has 
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not been found. This observation has been linked to the integrability of these 
models which restricts a full relaxation due to additional conserved integrals 
of motion. Nonetheless, numerical studies report that breaking integrability of 
strongly correlated fermions does not lead to relaxation [3T1I32] . and that the 
ID Bose-Hubbard model, which is commonly assumed to be nonintegrable, 
does not equilibrate for certain sets of initial conditions [33|3l] . 

This revitalized an ongoing debate about the long-term equilibration behavior 
of many-body quantum systems. It mimicked a similar discussion in nonlinear 
classical mechanics |35j which followed a seminal paper by Fermi, Pasta and 
Ulam [SH] • There the question was raised whether an excited system of coupled 
anharmonic oscillators would finally thermalize, i.e. approach a steady state 
thermal distribution of the initial excitation energy onto all oscillator modes. 
To the surprise of the original authors in 1955 this has not been confirmed by 
their numerical calculations. 

In the classical case, thermalization is regarded as a consequence of nonlin- 
earities in the equations of motion which generate chaotic behavior. Then, in 
general, the orbit of the classical motion samples completely the hypersur- 
face of the phase space that corresponds to all configurations which respect 
the conservation of certain extensive conserved quantities; then the system is 
called ergodicQ and thermalization is expected. Although the FPU problem 
is still debated [37], the predominant opinion today is that Arnol'd diffusion 
will finally lead to thermalization of the system. 

The time evolution of a quantum system is governed by the linear Schrodinger 
equation; it can always be represented by a strictly unitary operation U(t,t ) 
which is generated by the Hamiltonian. The later acts on the Hilbert space, 
such that an initial state is mapped onto a time-parametrized orbit of states 
constrained by constants of the motion like particle number, energy, etc. The 
time evolution of states is reflected in the evolution of expectation values of ob- 
servables; those correspond to measurable quantities. Note, for instance, that 
although the long time average of the dynamics of any quantum state equals 
its projection onto the zero energy (ground) state(sP] this does not indicate 



A classical system is called ergodic if two averages of the phase space coordinates 
coincide: a) The statistical (thermal) average over an ensemble of phase space con- 
figurations constrained to a hypersurface in phase space by constants of the motion 
(e.g. energy); and b) the time average of their dynamical solutions subjected to the 
same constants of the motion. This coincidence of a statistical (i.e. probabilistic) 
description of a physical system and a dynamical one (governed by deterministic 
laws of motion) is assumed by Boltzmann's hypothesis. Nonergodic cases can be 
easily found, e.g. systems with closed orbit solutions. It is an important result of 
nonlinear classical mechanics that nonlinearity in the equations of motions alone 
does not imply ergodicity. 

2 This can be easily seen by applying the time evolution operator in an eigenstate 
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a long-time relaxation towards (one of) its equilibrium ground state(s). 

Moreover, the strict linearity of quantum mechanical time evolution implies a 
fundamental disagreement between a quantum mechanical description based 
on the Schrodinger equation and a quantum statistical one by means of a 
statistical operator p; it is rooted in the additional coherence properties of a 
quantum system. Note that we only consider strictly closed quantum systems 
for which no tracing out of degrees of freedom related to an environment 
takes place. While a thermal state is constructed as an incoherent mixture of 
weighted quantum states with Tr[p^-] < 1, a quantum system, once initialized 
m a pure state defined by 1 = Tr[p 2 P ] = TT[U j {t,t )ppU{t,t )} remains pure at 
any later point in time because of the cyclic property of the tracd^J Due to 
the linearity of the trace, the same holds for any time average. Hence quantum 
systems are never ergodic in a classical sense. 

This limitation does not apply for the expectation value of a particular ob- 
servable. Nonetheless, the equilibration of expectation values in an integrable 
system is restricted for a different reason: Additional conserved integrals of mo- 
tion arise from an exact integration of the equations of motion which prevent 
a wipe-out of the influence of the initial conditions such that thermalization 
with respect to a conventional Gibbs ensemble should not be expected. In 
analogy to equilibrium approaches [3S], it was suggested that the long-time 
behavior can be reproduced by a statistical description based on a generalized 
Gibbs ensemble [26J . Many of the mentioned results [27|26|30j explicitly agree 
with this approach and its prerequisites and limitations have been discussed 
[39f24"] . A different notion of local relaxation grounds the examinations of fi- 
nite subsystems |4U|4*Tj which may exhibit thermal signatures even if the full 
system has not relaxed. 

Yet for closed nonintegrable many-body systems the fundamental questions 
became more obvious: Which observables exhibit thermal long-time behavior, 
and for what reason? How does the large number of degrees of freedom present 
in a many-body system make up for the unitarity of time evolution such that 
a thermal long-time limit can appear? Some earlier works which addressed 
these questions by introducing an eigenstate thermalization hypothesis [12115] 
found new attention recently [SJI5]. That hypothesis assumes that the ex- 
pectation value of a one-particle observable A in an energy eigenstate of the 
Hamiltonian H \M) = €m \M) equals the thermal average of the correspond- 
ing statistical quantity at the mean energy p = eu'- (M\A\M) = {A)(e M ). 
Note that this is an inherently time independent statement; it assumes that 
each single eigenvector of the Hamiltonian incorporates statistical signatures 



representation of the time-evolved state. 

3 Note that the cyclic argument is not applicable for expectation values of observ- 
ables (O) = tr[U^(t,t )ppU(t,t )O}. 
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and is sufficient to describe thermal behavior. A possible initial nonstatistical 
behavior appears as the result of a coherent superposition of eigenstates which 
dephases with time. 



1.3 Examination of quenched quantum systems 

In the following we will discuss a translationally invariant closed quantum 
system from first principle. The conservation of static quantities like energy, 
momentum, particle number constrain its dynamics; for any initial state it is 
deterministically determined at any time in the past or in the future by the 
Schrodinger equation. However, at zero time and by external influence, the 
Hamiltonian is changed abruptly. Additionally to a time independent (nonin- 
teracting) part Hq a two-particle interaction term Hi nt is switched on instantly 
such that 



0(t) is the Heaviside step function. To avoid trivialities, we assume that 
[H ,H int ] 7^ 0. This implies that the eigenbasis of the Hamiltonian changes 
at zero time: From the eigenvectors of H , {\m)} m , to the eigenbasis {\M)} M 
of the Hamiltonian H. For convenience, we will assume that nondegenerate 
perturbation theory can be applied; then a weak interaction does not change 
the noninteracting eigenenergies dramatically such that a relation can be es- 
tablished between the corresponding energies e m m €m for m = M. 

In general, a quantum quench implies that the system is open for the in- 
take of energy at the quenching time t = and for a redefinition of its 
Hamiltonian ground state. Yet it never means a sudden reset of its quan- 
tum state. Therefore, the quench initializes the interacting system described 
by the Hamiltonian H(t > 0) in the ground state of the noninteracting system 
\Q ) = \m = 0) which typically represents an excited state of the interacting 
system. Its properties can be studied by following its nontrivial successive dy- 
namics which is, again, generated by a time independent Hamiltonian. Hence, 
effectively, we study an initial value problem for a time independent Hamilto- 
nian H = Hq + -f^int- This observation determines our technical approach. 

Secondly, we choose a suitable observable to study the resulting dynamics. 
For translationally invariant systems the momentum mode number operator 
Mk is a convenient choice. It commutes with the noninteracting Hamiltonian 
[Af, Hq] = such that a common eigenbasis of both operators exists. This 
justifies the application of perturbation theory to study the dynamics of this 
observable^] As the momentum mode number operator is a one-particle ob- 




4 For this special situation the evaluation of the matrix element in the right hand 



H{t) = H + Q(t)Hi 



hit 



(1) 
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servable, it only exposes limited information on the interacting quantum sys- 
tem. This makes thermalization of its expectation value a plausible scenario. 



1-4 Solution of the Heisenberg equations of motion for an observable 



As a final part of the introduction, we construct an appropriate method to 
study analytically the effects of a quantum quench on observables. 

Working in a Heisenberg picture, states remain time independent while the 
operators carry the total time dependence. The dynamics of observables O 
which do not explicitly depend on time is described by the Heisenberg equation 
of motion 



m^^miH]. ( 2 ) 

This equation holds as an operator identity independent of any particular 
choice of a basis representation. Its solution, however, is best performed in an 
exact eigenbasis of the Hamiltonian where the dynamics of different energy 
scales is decoupled and can be treated separately. There time evolution leads 
to a pure dephasing of the Hamiltonian eigenmodes. 

In the following, we make this observation useful for approximating the time 
evolution of arbitrary observables in a more sophisticated application of per- 
turbation theory. Applying time dependent perturbation theory directly to an 
observable typically results in the generation of secular terms which are both 
proportional to time and the perturbative expansion parameter. They even- 
tually spoil the validity of time dependent perturbation theory even for small 
expansion parameters. Hence the aim is to separate a perturbative treatment 
of interaction effects from time evolution. This can be achieved by mapping the 
observable into an eigenbasis of the Hamiltonian where an exact evaluation of 
time evolution can be achieved. The corresponding transformation is the same 
which diagonalizes the Hamiltonian. Yet for many problems it is not known 
exactly. Implementing it by means of a perturbative expansion generates an 
approximate representation of the observable in an approximate Hamiltonian 
eigenbasis. Then secular terms arise again but, fortunately, in higher orders 
of the expansion only, such that the validity of time dependent perturbation 
theory has been improved. 

On the other hand, we have understood the quench scenario as an initial value 
problem and a competing requirement follows from the observance of nonequi- 



side of (??) is simplified; it could be approached, essentially, with a resummation 
of perturbation theory. Since for a fermionic system the eigenvalues of the number 
operator are bounded between zero and one, the result will depend on the filling. 
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Fig. 1. The Heisenberg equation of motion for an observable O is solved by trans- 
forming to the B = oo eigenbasis of the interacting Hamiltonian H (forward trans- 
formation), where the time evolution can be computed easily. Time evolution in- 
troduces phase shifts, and therefore the form of the observable in the initial basis 
B = (after a backward transformation) changes as a function of time. 

librium initial conditions. Note that the later are denned in the eigenbasis of 
the noninteracting Hamiltonian Hq. Particularly in the context of a many- 
body initial state, their mapping into a different basis representation may be 
complicated such that their evaluation is only convenient in their initial one. 

This motivates an approach which combines three unitary transformations: 
one which diagonalizes the interacting Hamiltonian approximately, called the 
forward transformation, a second one which represents a solution of ^ in 
the interacting eigenbasis and a final backward transformation (which reverts 
the diagonalization) into the noninteracting eigenbasis. All transformations 
are applied sequentially to the observable and implemented in perturbation 
theory. Note that the inverse transformation does not reproduce the original 
structure of the observable since time evolution has dephased the different 
contributions present in its forward-transformed representation. This scheme 
has been described in [13?] and is illustrated by fig. [I] It is motivated by 
canonical perturbation theory in classical mechanics. There the time reliability 
of a perturbative expansion can be greatly improved if it is performed after a 
suitable canonical transformation has been applied |48j . 



1.5 Outline 



In the following part of this paper we will apply the above approach first to 
systems with a discrete energy spectrum. As an introductory example we dis- 
cuss the exactly solvable model of a suddenly squeezed harmonic oscillator 
to illustrate technical aspects and the role of perturbative arguments. Then 
we will formulate a more general statement on the relation between equilib- 
rium and nonequilibrium expectation values for certain observables in discrete 
systems. We give two proofs which highlight different aspects of the character- 
istic nonequilibrium physics: The first proof stresses the role of perturbative 
arguments and restrictions imposed on the class of discussed observables by 
focussing on the overlap of eigenstates of a noninteracting Hamiltonian and 
its weakly perturbed counterpart. The second proof mirrors the operator ap- 
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proach depicted in fig. [T] It makes the drop-out of of transient and oscillatory 
behavior under time averaging more explicit. 

In a second part, we give details on the quench of a Fermi liquid as it has been 
described in pp. This constitutes a many-particle problem with a continuous 
spectrum of eigenenergies and requires more elaborate diagonalization tech- 
niques. Therefore we introduce the flow equation method following Wegner 
|19] which is an established approach towards the approximate diagonaliza- 
tion of many-body Hamiltonians and include it into the forward-backward 
transformation scheme. Although it is a nonperturbative technique we will 
only use it in an approximate form to set up perturbation theory. This is 
sufficient to observe the first phase of the dynamics of the quenched Fermi liq- 
uid. Studying the momentum distribution around the Fermi surface mirrors 
one-particle nonequilibrium physics. In particular, we observe a characteris- 
tic nonequilibrium value for the discontinuity of the momentum distribution 
at the Fermi surface which indicates the size of the quasiparticle residue. In 
equilibrium, correlation effects lead to a reduction of its value which is one 
only in the case of the noninteracting Fermi gas. In nonequilibrium, this re- 
duction is twice as large. This factor of two corresponds to the analogous 
perturbative results for the one-particle squeezed oscillator and in discrete 
systems. In a Fermi liquid, however, its appearance only indicates a transient 
nonequilibrium state. This transient behavior exhibits prethermalization [77] 
since contrary to the nonequilibrium momentum distribution the kinetic and 
the interaction energy have already relaxed to their final long time values. 
The relaxation of the momentum distribution gives rise to a second phase of 
the dynamics of a quenched Fermi liquid. It will be obtained from a quantum 
Boltzmann equation which describes the effective evolution of the momentum 
distribution from the nonequilibrium transient state onwards and leads to the 
prediction of its thermalization. 



2 Squeezing a one-particle harmonic oscillator 

The squeezed oscillator is a well-studied one-particle model system which 
found appreciation in many branches of physics. For two decades researchers 
have discussed squeezed states of the electromagnetic field which are interest- 
ing because of their characteristic reduced fluctuations in one field quadrature 
as compared to coherent states [oT)fo"T] . This suppression of quantum fluc- 
tuations in one variable out of a set of non-commuting variables below the 
threshold obtained for a state of symmetrically distributed minimal uncer- 
tainty, i.e. a coherent state, has motivated the naming: In this parameter the 
phase space portrait of the squeezed state shows sharp details and appears 
'squeezed' when compared to that one of a coherent state while fluctuations 
are inevitably increased in the others. The physical relevance of squeezed states 
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in optics is grounded on the fact that some interesting phenomena like gravita- 
tional waves generated in astronomical events are characterized by oscillations 
with amplitudes close to or below the width of the ground state wave function 
of an optical light mode as it is required by Heisenberg's uncertainty principle. 
Squeezed states, however, may provide improved signal-to-noise ratios beyond 
this quantum limit of coherent light and simplify the interferometric detection 
of the weak signatures of gravitational waves [53J. 

It has been shown [51] that squeezed states cannot be generated adiabatically 
from the ground state of a quantum mechanical oscillator but that sudden 
changes have to be applied to its parameters, e.g. its frequency or spring 
constant. Therefore, squeezed state represent an early example of what is 
now, in the context of a many-body system, called a quench of a quantum 
system. 

In many-body theory, the squeezing operation comes under the name of a 
Bogoliubov transformation. Recently, it was applied to study the behavior of 
a quenched Luttinger liquid in terms of bosonic degrees of freedom [27] . Here 
we use it to illustrate characteristic nonequilibrium behavior of one-particle 
models. Since the bosonic representation of a Luttinger liquid is momentum 
diagonal it equally serves as an example of effective one-particle behavior. 

2. 1 Hamiltonian 

On the level of the Hamiltonian, squeezing is inferred by an instantly applied 
change of the prefactor of the quadratic potential, namely the spring constant. 
We neglect a linear shift of the potential minimum and reduce squeezing to a 
sudden switch in the coupling constant g(t) = gO{t) of the quadratic particle 
non-conserving operators. With h = ujq = m = 1 

H = H + H int , H = o)a + ~ H int = g{t) ((a f ) 2 + a 2 ) (3) 

Representing the Hamiltonian in terms of space x = (a) + a)/ y/2 and momen- 
tum p = (a — a))/ y2 operators 

H=[l + g(t))x 2 +(^-g(t^p 2 (4) 

shows that it is strictly positive for \g\ < 1/2 and, thus, it is bounded from 
below. We only discuss the Hamiltonian for < g < 1/2. 

In the following we compare a perturbative analysis of this quench with an 
exact solution based on the exact diagonalization of H. In both cases we 
calculate the occupation, i.e. the expectation value of the number operator 
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N = crc both in the equilibrium ground state of the interacting Hamiltonian 
H and as a long-time limit of the dynamics of an initial state. We will show 
a remarkable relation between the equilibrium result and the nonequilibrium 
result which are described by the same functional dependence but variant 
prefactors. These nontrivial prefactors will later play a key role and can already 
be appreciated in this simple system. 

2.2 Perturbative study of squeezing 

Let us first assume that the coupling g is a small parameter and that nonde- 
generate perturbation theory can be applied to the interacting Hamiltonian. 
We formulate the perturbative approach in an operator language which cor- 
responds to the formalism of the flow equation method. 

2.2.1 Definition of the diagonalizing transformation 

The first step is to implement a discrete unitary transformation which diago- 
nalizes the Hamiltonian to leading order in g. We represent the unitary trans- 
formation U(<p) = e~ w by its generator rj, which is an anti-Hermitian operator, 
and by a scalar angle variable (p. Then the action of the transformation onto 
the Hamiltonian can be expanded according to the Baker-Hausdorff-Cambell 
relation as 

H = U^HU = e^He-w « H +H- mt +(fa H ] + [v, H int ]) ^+-[77, [ V , H ]] V 2 +. . . 

Demanding that to leading order the interaction term should vanish leads to 
an implicit definition of the generator 

ip[rj,H ] = -H iMt . (5) 

Note that it implies rj ~ 0{g). It can be easily checked that the canonical 
generator defined as the commutator of the noninteracting and the interacting 
part of the Hamiltonian fulfills this implicit definition ([5]) if an angle <p=\/A 
is chosen: 

V = [H ,H iat ] =2g((a r ) 2 -a 2 ) (6) 

2.2.2 Transformed Hamiltonian 

In a second step we consider the corrections beyond leading order in the (ap- 
proximately) diagonalized Hamiltonian which are, in general, second order in 
9- 
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H = H + [ Vl H mt }ip+ 2 [r ] ,[r ]j H ]]ip 2 = H -16 V g 2 H + 32 V 2 g 2 H + O(g 3 ) 

^ /A (l-2g 2 )H + O(g 3 ) 

Thus the transformed Hamiltonian is described by a renormalized frequency 
oj — (1 — 2g 2 )Q Due to the particular simplicity of squeezing a harmonic os- 
cillator the second order correction can be fully absorbed in a renormalization 
of parameters. 



2.2.3 Transformation of quantum mechanical observables 

Similarly to its action onto the Hamiltonian the unitary transformation implies 
a transformation of all quantum mechanical observables which constitutes the 
third step of a unitary diagonalization approach. 



6 = + [r,,0]^+^[ V , [ V ,0}] V 2 + ... 



(7) 



We make the transformation of creation and annihilation operators explicit 
up to second order in g and write 



l+9 2 /2 -g 
-9 1 + W2 



== T{g) 




The three steps (2.2.1 - 2.2.3) establish a diagonal representation and are, 



altogether, referred to as the forward transformation. We note that the trans- 
formation of the observables can be easily inverted. Up to second order in g the 
inverse of T, called the backward transformation, is given as T~ l (g) = T(—g). 



2.2.4 Spin-off: The equilibrium occupation 

We interrupt our calculation of the nonequilibrium occupation for a short 
detour in order to evaluate the equilibrium one. To be more specific, our in- 
terest is in the occupation of the interacting ground state with 'physical' par- 
ticles (i.e. particles which are defined by the eigenmodes of the interaction-free 
Hamiltonian Ho). Denoting the interacting ground state by \Q) the equilibrium 
occupation reads iV E Q u = (0| a^a \Q). This is unitarily equivalent to the eval- 
uation of a transformed number operator with respect to the noninteracting 
ground state |O ) since iV EQU = (0| WUdaJUfiU |0) = (fi„|WoW|fi). For- 
tunately, the transformation which links both representations is the inverse of 

5 We point out that this equation constitutes the discrete analogue of a flow equa- 
tion for the Hamiltonian (cf. ??). 
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the forward transformation. Hence we evaluate with n a = (Qq\ a' a \ Qq) up to 
second order in g 



( 2 \ 2 

N EQU = g t 5 |^ ^ / 1 + 9_ \ fl t a ^ + g 2 fla t ^ 

^(l + 2( 7 2 K + ( ? 2 (9) 

We note that the occupation of the oscillator measured in terms of the original 
'particles' is increased. In the following we will compare this result with the 
nonequilibrium occupation obtained after sudden squeezing. 



2.2.5 Time evolution of transformed observables 

We resume the calculation for the nonequilibrium case. The forward trans- 



formation has already been completed in (2.2.1 - 2.2.3). In a fourth step the 
transformed observables are time evolved for all positive times with respect 
to the transformed Hamiltonian. This, effectively, accounts for the insertion 
of time dependent phase factors. 




e iHt~^ e -iHt 



e iHt~ e -iHt 



e~ iuJt a 



a e^(l + g 2 /2) 



_ e -iut g 



-e iuJt g 



:i+W2) 



) 









2.2.6 Backward transformation 

Finally, we map back the time-evolved observables to the eigenbasis of the 
noninteracting Hamiltonian, completing the scheme depicted in fig. [T] Up to 
second order in g we obtain 




T-\g) 




e iut + 2igZ s in(W) -2ig(l + g 2 /2) sm(ut) 
v 2ig(l + g 2 /2) sm(ut) e~ iuJt - 2ig 2 sm(ut) 



\ 



This constitutes a consistent perturbative solution of the Heisenberg equations 
of motion for the operators a) and a. 
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2.2.7 Nonequilibrium occupation 

In a final step we compose the time dependent number operator from the time 
dependent creation and annihilation operator in an obvious way. Since time 
evolution is unitary, the time evolution of a product of operators is always the 
product of the time evolved operators which can be easily checked by inserting 
unity H = U(t,t )U Jl (t,tQ). Evaluating the expectation value of the number 
operator for the initial state \Qo) leads to the nonequilibrium occupation 

N NEQ (t) = (tt \a^(t)a(t) \n ) = n a + Ag 2 (2 sin 2 (ut))((n \a^ a \n ) + i) (10) 

The large time limit is obtained by time averaging which is defined for a time 
dependent variable A(t) as A :— ^ J T dtA{t). Then j\r NEQ = n a + Ag 2 n a + 2g 2 . 
Comparing with ([9]), we find with AN(t) := N(t) — n a 

lim AiV NE Q(t) = 2 AN EQU + Og 3 (11) 

t — >oo 

The factor of two between the equilibrium and the nonequilibrium occupation 
constitutes the main result of this calculation. It states that even in a long-time 
limit the nonequilibrium occupation does not approach the equilibrium one. 
The numerical value of two can be considered as a consequence of applying 
two transformations, the forward and the backward one, such that changes to 
the occupation due to interaction effects double. In the following we will find 
that, although the numerical value gets corrections in order g 3 , the mismatch 
of both occupations is retained for all orders of perturbation theory. 

2.3 Nonperturbative (Bogoliubov) treatment of squeezing 

In a second approach, we implement a Bogoliubov transformation which ex- 
actly diagonalizes the squeezing Hamiltonian. It can be found in many text- 
books, e.g. [5Tf52] . Our aim is to illustrate that the chosen perturbative ap- 
proach exhibits, up to numerical details, the correct nonequilibrium behavior 
of the system. 

The exact diagonalization of the squeezing Hamiltonian ^ can be constructed 
from the action of the (inverse) unitary squeezing operator 

where £ = re %e is an arbitrary complex number which will be specified later. 
Applying the squeezing operator to the ground state generates squeezed states 
in analogy with the displacement operator which maps the ground state onto 
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coherent states. Here we go the opposite way and apply its inverse to diago- 
nalize the squeezing Hamiltonian. 



2.3.1 Exact transformation of observables 

We directly start with writing down the action of S(£) onto the creation oper- 
ator and the annihilation operator. In condensed matter theory it is commonly 
known as a Bogoliubov transformation used to treat interactions quadratic in 
creation or annihilation operators. 



at 



s(0 



( 
v 



cosh(r) — e %e sinh(r 
-e sinh(r) cosh(r) 




T F {i) 



a' 
a 



(12) 



Note that T F (£) is not a unitary matrix despite the fact that det(T F (£)) = 1. 



2.3.2 Exact Hamiltonian diagonalization 

Inserting this transformation into the interacting Hamiltonian ^ results in a 
sum of four terms: 



H 



f^) 2 x 



e ld cosh(r) sinh(r) + g cosh 2 (r) + ge 2td sinh 2 (?/ 

cosh 2 (r) + sinh 2 (r) + 4gcos(#) cosh(r) sinh(r) 
11 x [2 cos(#) cosh(r) sinh(r)] 



a) a x 



+ h.c. + 
+ 



To achieve a diagonal Hamiltonian we demand that the terms quadratic in 
a) and a should vanish. This fixes the free parameter £ := re 10 = r. For 
small interactions \g\ < 1/2 real solutions with 8 = can be found. With 
sinh(2r) = 2sinh(r) cosh(r), cosh(2r) = sinh (r) + cosh (r) the real parameter 
r can be linked to the interaction 

tanh(2r) = -2g 

For small values of g <C 1/2 the expansion arctanh(x) ~ x for i< 1 implies 
that r w —g. Then the nonperturbatiye Bogoliubov transformation coincides 

VP=l/4' 

cosh(2r) + 

25fsinh(2r) compared to the original frequency ujq = 1 in Q3p . For all values of 
g < 1/2 the renormalized frequency is positive and the Hamiltonian is bounded 
from below. Its dependence on g is plotted in fig. [2j In the limit of small g we 
find for the renormalized frequency its perturbative value to = (1 — 2g 2 ). 



with the perturbative approach in |2.2 ), with S(£(g) 
The diagonal Hamiltonian shows a renormalized frequency u 
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2.3.3 Exact equilibrium occupation 

Again, we first calculate the expectation value of the equilibrium number 
operator, using cosh (r) — sinh (r) = 1 



N EQU = ^| a ] a |^ = ( 0q | -tg |^ 

= cosh 2 (r) (f2 | a) a |Q ) + sinh 2 (r) (fi | aa) \Q Q ) 

= n a + 2sinh 2 (r) (n a + -) (13) 



Again, the perturbative limit for small g agrees with ([9]). 



2.3.4 Exact nonequilibrium occupation 



For the nonequilibrium occupation we solve the Heisenberg equations of mo- 
tions for the creation and annihilation operators in the (now exact) eigenbasis 
of the Hamiltonian. The forward transformation of these operators is given by 
(12). Again, we complete the scheme in fig. [I] and compute the time evolution 
of the transformed operators with respect to the diagonalized Hamiltonian, i.e. 
with respect to the renormalized frequency u. The final backward transforma- 
tion is given by T B (r) = T F (—r). These three steps can be easily denoted as 
subsequent matrix multiplications: 



=T B (r) 



T b (r) 



e 



itut I 



cosh ( r ) — e sinh ( r 



i sin(co>t)e sinh(2r 



^ — i sm(iot)e ie sinh(2r) — (V^cosh 2 



-iuit 



sinh (r 



cos(ut) 



1 
-1 



— i sin(u;i) 



cosh(2r) — e sinh(2r 
1 e ld sinh(2r) cosh(2r) 




Composing the number operator reads 



N NEQ {t) = (Q \a\t)a{t) |0 Q ) 



(Q \ [(e-^ cosh 2 (r) - e iuJt sinh 2 (r)) 0) + i sm(tot)e- ie sinh(2r)a 



-ism(ut)e ie smh(2r)a j - (e iw *cosh 2 (r) 



-iuit 



sinh (r)) a 
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Fig. 2. This plot illustrates the nonperturbative ratio m{r{g)) 



sinh 2 (2r(g)) 



in (15) 



2sinh J (r(g)) 

between the equilibrium and nonequilibrium correction to the noninteracting occu- 
pation due to interaction effects. The real solution for the diagonalization transfor- 
mation £(c/) = r{g) is valid only for g < 1/2. In the limit of small interaction the 
factor m = 2 is exact. Additionally the renormalized frequency u>(g) of the diagonal 
Hamiltonian is shown. 

Only the particle number conserving terms ~ a^a, aa) contribute and we arrive 
at the nonequilibrium occupation 



N NEQ {t) =n a + (2sin 2 (u;t)) sinh 2 (2r) 



n a + 



Again, the long time limit is taken as a time average. This implies that the 
renormalization of the frequency does not affect the occupation at late times 
and may be neglected. 



lim N NE Q(t) =n a + sinh 2 (2r) (n a + 1 



(14) 



2.3.5 Nonperturbative relation between the equilibrium occupation and the 
nonequilibrium occupation 



Comparing (14) with (13) one observes that for the squeezing Hamiltonian 



the relation between the equilibrium and nonequilibrium occupation is given 
by 



m(r(g)) := 



AN neq ^ S inh 2 (2r(#)) 
AATEQU ~ 2sinh 2 (r(^)) 



(15) 



The precise numeric value of the ratio depends via r on the coupling strength 
g and is plotted in fig. [2j Expanding sinh(r) « r + r 3 /3! + 0(r 5 ) we confirm 
that in the perturbative limit this relation approaches a factor of two. 
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3 More general statements on quenched one-particle systems 



As the discussion of the squeezed oscillator has shown a perturbative approach 
to quenched one-particle systems captures important signatures of the behav- 
ior of their nonequilibrium occupation. In the following we will show that it 
can be generalized to a large class of weakly interacting one-particle systems 
and to more general observables. 

Prerequisites: Let us consider the Hamiltonian of a quantum system with a dis- 
crete energy spectrum and small interaction g which models a weak quantum 
quench. 

H = H + gQ(t)H int (16) 

Its interacting ground state is denoted by |0). We assume that nondegenerate 
perturbation theory with respect to the noninteracting ground state |O ) is 
applicable and that Hq and H mt do not commute. 

Moreover, we assume a quantum mechanical observable O which does not 
depend explicitly on time and obeys the following relations: 



£>|O )=0 (17) 
[O,H ] = Q (18) 

For times t > its time evolution 0(t) is generated by the interacting Hamil- 
tonian Hq + gH mt . Then the following theorem holds: 

Theorem: In second order perturbation theory the long-time limit of the time- 
averaged expectation value of the time evolved observable in the initial state 
equals two times the equilibrium expectation value of the observable in the 
interacting ground state. 

O : = lim i [ T dt (O | 0(t) |O )=2 (01 O 10) + 0(g 3 ) (19) 

T^oo T JO 



Proof 1: We will give two proofs of this theorem. The first one is intended to 
motivate the physical origins of the prerequisites by relating it to the more 
conventional picture of overlapping eigenstates. This allows to conclude on its 
general relevance. 

Firstly, we introduce eigenbasis representations for the noninteracting Hamil- 
tonian H {\m)\m G -/V }[][] the interacting Hamiltonian {|M)|M 6 A^pjf 7 ] 

6 We label eigenstates of the noninteracting Hamiltonian by lower case variables. 
We label eigenstates of the interacting Hamiltonian by upper case variables. 
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and the observable { \ E N } with the eigenvalues e° m = €m + 0{g 2 ) (for 



m = M) and Oj, respectively. The requirement (18) implies the existence of a 
common eigenbasis of the observable O and the noninteracting Hamiltonian 
Hq such that we can assume pairwise coinciding eigenvectors \m) = For 
clarity, however, we will keep a separate notation. The equilibrium ground 
state expectation value is rewritten by inserting unity. 

(Q\0\Q) = J20 J \(j\Q)\ 2 (20) 

3 

An analogous evaluation of the time dependent expectation value by inserting 
unities, extracting time dependent phase factors and taking their time average 
leads to 



(Q \e~ iHt Oe iHt \n ) = 

= lim i f T dt £ (fi | e~ im \M) (M \j) (j\ O \f) (f \M) (M'\ e iHt \Q ) 

T-,00 T Jo MM'jj' 

= lim i [ T dt ^ Oi/e^—)* (O | M) {M \j) (j\f) (f \M) (M'\n ) 

T-oo T Jo MM'jj' 

= J2O J \(M\n )\ 2 \(j\M)\ 2 (21) 

Mj 

Up to a relative phase, the interacting eigenstates \M) are invariant under 
time evolution. Therefore, overlap matrix elements are discussed with respect 
to these states. 

The set of matrix elements {(M \ flo)}M describe a decomposition of the initial 
state in terms of Hamiltonian eigenstates. This is a statement about the partic- 
ular initial conditions of the quench problem. Since we discuss a quench from 
the noninteracting Hamiltonian this is a decomposition of the noninteracting 
ground state in terms of interacting eigenstates. 

The second set of matrix elements {(j \M)} M j encapsulates the overlap be- 
tween the eigenbasis of the observable and the eigenbasis of the Hamiltonian. 



Since (18) holds we can work in the common eigenbasis of the observable and 
the noninteracting Hamiltonian. Then the overlap between the eigenbasis of 
the interacting and of the noninteracting Hamiltonian is discussed. In both 
cases the matrix elements can be evaluated by applying perturbation theory 
to the Hamiltonian H = Ho + H[ nt , treating H\ nt as a small perturbation. We 
make this explicit to leading order: 

for M = m 
(mfnM 2 for M ^ m 



\(m \M}\ 2 = 



As (17) implies O = the direct overlap between the interacting and the 
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noninteracting ground state does not contribute to the sums in both (20) and 



(21 ); hence they are at least second order in g. We compare the right hand side 
of both equations for any fixed value of j. In the nonequilibrium case, second 
order contributions require a resonance condition for the involved quantum 
numbers, M = j or M = 0. 



|(M|fi )| 2 x|(j|M)| 2P = T 



\(J \n }\ 2 x 1 for M = j =: J 
lx\(j\n)\ 2 forM = 



Because of the symmetry \(J |fio)| 2 = \(j |^)| 2 i n leading order perturbation 
theory, both contribute equally \(j |fi)| 2 to the sum over M. Then in second 
order perturbation theory holds 



(n \e- iHt Oe im \Q ) =2j20 j \(j |fi)| 2 = 2 (Q\ O 



and the theorem is proven. 

Proof 2: The second proof of the theorem is constructed in analogy to the 
scheme presented in fig. [T] and aims at a clearer understanding of the applied 
method and its particular merits. 

(1) Definition of a unitary transformation. We define a single unitary trans- 
formation U\ = e Va by its anti-hermitian generator t] s = —T]\, demanding that 
its application to the Hamiltonian disposes the interaction part of the Hamil- 
tonian to first order of g. Expanding its unitary action onto the Hamiltonian 



H = e^He-** = H +H int + [r) s , H )+[r, s , H iQt } + - [r, s , [r, s , H + H iQt ]] + 0(g 3 ) 

(22) 



O(g) 



hence allows to read off an implicit definition of r) s by 

[rj 8 ,H ] = - H int 



which justifies the assumption r] s ~ O(g) in (22). Then the transformed Hamil- 



tonian equals the free Hamiltonian up to second order corrections. 

(2) Computation of the interacting ground state expectation value. In the fol- 
lowing we exploit a formal coincidence which holds for all systems with a 
nondegenerate single ground state: For any such Hamiltonian, the diagonal 
representation of the interacting ground state in terms of the diagonal degrees 
of freedom can be formally identified with the ground state of the noninteract- 
ing Hamiltonian; thus we can relate them by li\ |fi) = |O ) or, to leading order, 



19 



by \ = (1 —77s) |fio) + 0(g 2 ). As every Hamiltonian can be diagonalized, this 



does not pose any further restrictions. Hence with (17) 



(n\ o \q) = (n \u$ou s \n ) = -(n \ Vs Or] s \n ) + o( g 3 ) (23) 

The simple diagonal representation of the interacting ground state motivates 
the application of operator-based transformation schemes like the flow equa- 
tion method in equilibrium since correlation effects are, formally, fully trans- 
ferred from the description of an interacting ground state to the the particular 
form of transformed observables. Thus one can avoid to discuss the full com- 
plexity of the interacting ground state and restrict to those correlation effects 
which become actually relevant for a particular observable. The transforma- 
tion U\OU s can be performed in the most convenient way. 

(3) Real-time dynamics of the observable after the quench. For the evalua- 
tion of the nonequilibrium expectation value O we start with the sequential 
application of three unitary transformations. Firstly, at time t — the ob- 
servable is represented approximately in the energy- diagonal eigenbasis of the 
Hamiltonian. 

6(0) = O(0) + [77,, O(0)] + l - [77,, [ Vs , O(0)}} + 0(g 3 ) (24) 

Now we apply unitary time evolution to the transformed observable with re- 
spect to H = H + 0(g 2 ). This is time-dependent perturbation theory to first 
order. 

Q(t) = e- iHot O(0)e iHot (25) 



We insert (24) into ([25]) and attribute the time dependence to the generator 



Vs — * Vs(t) — e tH ° t r] s (0)e' lHot . This is possible because of (18) and ensures that 



(17) holds for all times. Finally, the backward transformation t)b = —rj s (0) is 



applied. 

0(t) = &(t) - [77,(0), 6(t)} + ~ [77,(0), [77,(0), 0(t)}} + 0(g 3 ) 
= O(0) + [ Vs (t), 0} + \ [ Vs (t), [ Vs (t), O}} - [77,(0), O} 

- [77,(0), [ Vs (t), O}} + \ [17.(0), hit), O}} + 0(g 



3\ 



We evaluate the expectation value of 0(t) in the initial state |f2 ). Due to (17) 
many contributions vanish. 

(Q \ 0(t) \Q ) = (2 (-1) 77 s (t)0r/ s (t) + r/ s (0)Or/ s (t) 



2 



+77 s (t)0r/ s (O) - 2 ( I ) Vs (0)O Vs (0) ) (26) 



n 
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Inserting unity 1 = J2 m \ m ) ( m \ i n terms of eigenstates of the noninteracting 



Hamiltonian H shows that the second and the third term in (26) dephase 
and do not contribute to the long time average: 



(fi | r) s (0)Or) s (t) \Q ) = £ W rj.(0)O \m) (m\ e~ lHot r] s (0)e lHot |fi > (27) 

m 

= EO m e l(£ ™- £0)1 (n \ 77.(0) \m) (m\ r/ fl (0) |fl„) (28) 



For m = equation (17) implies Oq = 0. As we have assumed a nondegenerate 
Hamiltonian H we obtain 

(n \ Vs (o)Or) s (t) \n ) = (fi | vs(t)o Vs (o) \n ) = o 



On the other hand, making use of (18) 



(fi | e-^^e^^e-^e^ * |O ) = (Q | ^(0)0^(0) |O ) 



Consequently, we arrive at C = -2 (fi | r/ s (0)0^(0) |n >- With (23|) the the- 
orem is proven. 

The second proof explains the factor of two as the accumulation of equal 
second order corrections both from the forward and from the backward trans- 



formation. The drop-out of transient or oscillatory behavior in (26 ) due to time 
averaging is more explicit. This depicts the major merit of the transformation 
scheme: Fundamental correlation-induced effects - as it is, for example, the 
difference between the interacting and the noninteracting ground state- enter 
a perturbative study of time evolution performed in an energy-diagonal repre- 
sentation already as time-independent offsets. That their influence is stronger 
in nonequilibrium than in equilibrium can be seen directly. 

Corollary: In many systems the noninteracting part Hq of an interacting 
Hamiltonian H represents the kinetic energy for which the following relation 

holdsfj] £NEQ, KIN _ 2£'EQU,KIN_ 

Proof: Define £ NEQ > KIN (t) := (fi | H Q {t) |fi ) and apply the theorem for O = 
H . 



Discussion: 

The above theorem explains the factor 2 in the ratio between nonequilibrium 
and equilibrium expectation values as a rather general observation in systems 



8 The increase in the kinetic energy beyond its value for the interacting ground 
state, however, does not necessarily indicate its thermal distribution (i.e. heating 
effects). 
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with discrete energy spectra. In the following we will see how this factor 2 
appears for an interaction quench in a Fermi liquid with continuous spectrum 
and what role it plays for the nonequilibrium dynamics. Notice that recently 
similar observations have been made for a quenched Kondo impurity in the 
ferromagnetic regime [55]. There the interaction of the impurity spin with a 
band of metal electrons is switched on suddenly in time and the subsequent 
spin dynamics is calculated. The nonequilibrium dynamics of the magnetiza- 
tion has been studied by both analytical (flow equations) and exact numerical 
(time dependent NRG) methods [55]. The respective results agree very well 
on all time scales and again show the above factor 2 when comparing to the 
equilibrium value. 



4 Real-time evolution of a quenched Fermi liquid 

In the last section we have compared the equilibrium and nonequilibrium be- 
havior of a single-particle anharmonic oscillator. In the following we present 
an analogous analysis for a many-body noninteracting Fermi gas for which 
we compare a sudden switching-on of a two-particle interaction with an adia- 
batic evolution into an interacting ground state. In more than one dimension 
the later corresponds to Landau's theory of a Fermi liquid, which, since its 
introduction in the late 1950s [55], became a benchmarking effective descrip- 
tion for the study of many (normal) interacting Fermi systems [57.58J. Its 
main prerequisite (at least from the point of view of this work), is the adi- 
abatic connection between the noninteracting Fermi gas and the interacting 
Fermi liquid. It means that there is a continuous evolution from the low en- 
ergy states of the gas to those of the liquid as the interaction is increased. 
This continuous link allows to formulate the intuitively appealing concept of 
Landau quasiparticles since particle properties carry over from the noninter- 
acting physical fermions to the interacting degrees of freedom: Around the 
Fermi surface, the interacting degrees of freedom, then called quasiparticles, 
differ from the noninteracting ones only by modified parameters, e.g. an effec- 
tive mass. In terms of physical fermions Landau quasiparticles are composite 
many-particle objects although they are not true eigenstates of the system; 
this implies both a residual interaction among them and their finite lifetime 
which roughly measures the departure from interacting eigenstates. As the 
lifetime diverges right at the Fermi energy Landau's approach becomes ex- 
act there. Away from the Fermi energy two incompatible time scales compete 
with each other, limiting the applicability of Landau's approach to low energy 
excitations: The adiabatic requirement of long ramp-up times and the limited 
lifetime of the achieved quasiparticle picture. In this work we will discuss the 
opposite limit of the adiabatic increase of the interaction strength by applying 
a sudden quench. We will find that for weak quenches within the Fermi liquid 
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regime a quasiparticle picture can be retained. 
4-1 Hubbard model 

We discuss a generic Fermi liquid by referring to the Fermi liquid phase of the 
one-band Hubbard model [I] ; its equilibrium properties have been extensively 
studied using a great variety of different methods [59; 60f61U62f63f64j . It repre- 
sents a lattice model which we discuss in the thermodynamic limit, i.e. for an 
infinite number of lattice sites N. We assume translation invariance. Firstly, 
we do not specify a particular lattice geometry, work with a generic dispersion 
relation e(k) but assume that no nesting of the Fermi surface occurs. Each 
lattice site is capable of carrying up to two spin - 1/2 fermions. Due to the 
Pauli principle this implies a local state space of dimension fouip~[ Without 
interaction, all local states are degenerate in energy; hence their composition 
to the state space of a many-site lattice leads to a single energy band. We 
discuss the model for finite bandwidth D and in the regime of half filling, i.e. 
with -on average- one fermion per lattice site. The Hubbard Hamiltonian 



displays two competing physical processes on the lattice: The first term de- 
scribes cohereni P^] hopping processes between two neighboring lattice sites j 
and i with site-independent strength = t = 1 which we set equal to one for 
convenience. This defines our energy scale in the sequel. The second term de- 
picts a site-independent repulsive interaction U. It approximates the influence 
of a Coulomb repulsion of electrons dwelling on the same lattice site. Due to 
the Pauli principle, it is effective only between electrons of different spin and 
is proportional to the product of their spin densities. It reduces the mobility of 
fermions on the lattice as, intuitively speaking, 'hopping onto singly occupied 
sites becomes energetically less favorable'. Hence the Hubbard model depicts 
the competition between delocalizing and localizing effects in an interacting 
Fermi gas. At zero temperature a Fermi liquid phase exists for weak interac- 
tion strength in all dimensions larger than one. The chemical potential /x is set 
to half filling. In a momentum representation it merely accounts for a global 
energy offset of the kinetic energy of the fermions. The allowed momenta are 



Note the following configurations of the four-dimensional local state space: zero 
occupation, the spin up or spin-down configuration of single occupancy and the 
combination of two antiparallel spins (double occupation). 

10 An intuitive but not fully correct imagination depicts these hopping processes 
as an exchange of fermions between different lattice sites. Yet the coherent nature 
of the hopping is better captured by their delocalizing effect on a single fermionic 
wavefunction which spreads over more than a single lattice site. 




(29) 
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restricted due to the limited bandwidth /C = {k 6 R d : e{k) G [— D,D]}. For 
simplicity, Arabic numbers are used as generic momentum index labels. Note 
that the Hubbard interaction, while being local in real space, is non-diagonal 
in momentum space. 

H = Yl e k c Lcka + u Yl ct i'T c it4'|C2i ^Kf' ( 30 ) 

fceK, CT e{U} i'i2'2eK 



4-2 Quench of a Fermi liquid 

We implement a particular interaction quench of a Fermi liquid. It is modeled 



by substituting U U(t) = G(t)U in (30). 



4-2.1 Normal ordering and energy considerations 

We decompose the interaction term by applying a normal ordering procedure 
with respect to the ground state of the Fermi gas \fl ) = |FG), which equals 
the noninteracting Fermi gas (FG). Its momentum distribution is given by 
nk := (f^ol cjxfc |Oo) and does not depend on spin. Normal ordered operators 
are denoted between colons. 



t t t t r2' t rl' t rl' s-2' 

c i'T c it c 2'j. c 2| = :cj, T Ci|C2/jC 2 j: + n 2 >o 2 :c[,^Cif. + n v b l :c 2 ,^c 2 i: — n v n 2 'O x d 2 

Thus one particle properties hidden in two-particle scattering terms are ex- 
tracted and a clear separation of one- and two-particle features in the Hamil- 
tonian is achieved [68]. The one-particle scattering contributions correspond 
to a shift of the chemical potential of — U/2 and could be reallocated to the ki- 
netic energy. However, this additional energy is not dynamically relevant. We 
stress the point that the observed dynamics is solely caused by two-particle 
interaction effects. 

Then the sudden switch-on of a normal ordered two-particle interaction term 
does not lead to a change of the total energy of the system. Yet it lowers 
the ground state energy of the Hamiltoniar j 11 1 such that, at zero time, the 
noninteracting ground state of the initial Fermi gas, |fio) ; is promoted to a 

11 This consequence of an energy conserving quench may seem unphysical, but it 
highlights a simple fact: While adding a (Coulomb) charge to formerly uncharged 
particles leads to an additional one-particle potential energy which is first order in 
U and can be accounted for by a shift of the chemical potential, the corresponding 
two-particle repulsion (which causes interparticle correlations and lowers the ground 
state energy) is a second order effect; the dynamics of the quenched Hubbard model 
is only driven by the later. 
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highly excited state of an interacting model. The corresponding excitation 
energy is measured with respect to the ground state energy of the Hubbard 
Hamiltonian in equilibrium. It does not vanish in the thermodynamic limit as 
a single-particle excitation would do. This is because the correlation-induced 
reduction of the ground state energy becomes effective at every lattice site. 
A more detailed discussion of the involved energies will be resumed in section 

E3U 



4-2.2 Observables 

To characterize the dynamics of the quenched Fermi liquid we analyze the 
evolution of particular quantities, namely the total kinetic energy, the total 
interaction energy and the momentum distribution function. They are expec- 
tation values with respect to the initial state \ fl ) of the observables H , H int 
and the number operator for a fermionic quantum gas 



We work in a Heisenberg picture where the observables carry the time de- 
pendence. Notice that the distribution function exhibits the evolution of the 
occupation of one-particle momentum modes while the energies are mode- 
averaged quantities. As the kinetic energy can be easily calculated from the 
momentum distribution function; and as the total energy of a closed system 
is conserved, it is sufficient to explicitly calculate the momentum distribution. 

In equilibrium, the zero temperature momentum distribution of an interact- 
ing many-particle system of fermions is characterized by a discontinuity at 
the Fermi momentum kp. Its size, the so-called quasiparticle residue or quasi- 
particle weight < Z < 1, reflects the strength of interaction effects: For 
an interaction-free Fermi gas it acquires its maximum value one, while it de- 
creases with increasing interaction strength in a Fermi liquid. For the equilib- 
rium Hubbard model this signature has been studied numerically throughout 
and beyonc p^] the Fermi liquid phase [SofflT] . The behavior of the quasiparti- 
cle residue under nonequilibrium is the main focus of this work. It is, clearly, 
a zero temperature analysis as temperature smears out the discontinuity on 
its own energy scale. 



At a critical interaction strength, marking a quantum phase transition from the 
Fermi-liquid to a non-metallic (Mott insulator) phase, the quasiparticle residue van- 
ishes. 




c[,Cfc : k > kp 
1 — c\ck : k < k F 
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5 Perturbative analysis of a quenched Fermi liquid 



Similarly to the perturbative study of the squeezed harmonic oscillator we 
implement the time evolution of the number operator following the scheme of 
fig. [T] for a Fermi liquid. Again we aim at a perturbative analysis, expanding 
all results as a power series in the interaction strength. Contrary to the single 
mode oscillator, the Fermi gas in the thermodynamic limit represents a many- 
particle problem with an infinite number of degrees of freedom. This implies 
that many different energy scales contribute to the Hamiltonian. It is obvious 
that the quench generates occupation in many different excited eigenstates of 
the interacting Hamiltonian. Therefore we implement the diagonalizing trans- 
formation such that a controlled treatment of different energy scales is possible. 
This can be achieved by a flow equation transformation following Wegner |1H] . 
Due to the Pauli principle, a fermionic many-particle problem is characterized 
by the existence of a filled Fermi sea; the later restricts the phase space for 
fermionic scattering processes, in particular at zero temperature. This, effec- 
tively, reduces the strength of the two-particle Hubbard interaction and allows 
for the observation of a transient dynamics of an excited state. 



5.1 Flow equation transformation 



Since it has been independently introduced by Wegner jl9] and Glazek and 
Wilson [66.67] the flow equation method became an established tool in the 
analysis of equilibrium and nonequilibrium many-body systems and has been 
applied to a great variety of different systems. An extensive list of model 
systems and problems which have been tackled by the flow equation method 
can be found in [HE] and a comprehensive textbook review is available [69] . 
Quite recently, it has been successfully applied to nonequilibrium problems 
[70|7T|fT] . For the convenience of the reader we will give a short introduction 
here. 



5.1.1 Continuous sequence of infinitesimal transformations 



In section (2.2.1 ) a single unitary transformation was defined to diagonalize the 
Hamiltonian approximately. The flow equation method, however, decomposes 
the diagonalization of a many-body Hamiltonian into a continuous sequence 
of infinitesimal unitary transformations. It aims at an approximate diago- 
nalization of the Hamiltonian in energy space. Applied to the nondiagonal 
Hamiltonian, it imposes a continuous evolution of Hamiltonian parameters. 
This evolution is, in rough analogy to Wilson's interpretation of the renor- 
malization group, depicted as a flow of the Hamiltonian parameters. The flow 
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is parametrized by a scalar, nonnegative and monotonously growing flow pa- 
rameter B. 

Since only the initial Hamiltonian and, to a lesser degree, the final, approxi- 
mately energy- diagonal Hamiltonian are fixed boundary conditions there is a 
large degree of freedom how the continuous sequence of infinitesimal unitary 
transformations is actually constructed. It allows for the implementation of 
other desirable features like energy scale separation such that the flow pa- 
rameter B can be related to an energy scale Ab = l/y/B. But contrary to 
conventional renormalization schemes which distinguish between absolute en- 
ergies of high and low energy degrees of freedom the flow equation methods 
separates the treatment of large and small relative energy differences in the 
Hamiltonian. This means that those matrix elements of the Hamiltonian which 
describe deeply inelastic scattering processes are eliminated already at an early 
stage of the flow. Successively, those with lower energy differences are treated 
while elastic scattering processes ('energy-diagonal ones') remain unchanged. 
Hence, the flow equation method achieves a stable sequence of perturbation 
theory. As all energy scales of the Hamiltonian are retained this motivates the 
application of the flow equation method to nonequilibrium phenomena. 



5.1.2 Definition of the infinitesimal transformations 

Wegner showed that energy scale separation can be implemented by defining 
the canonical generator of infinitesimal transformations [49], representing a 
differential form of equation (|6]). 



r)(B) = [H (B),H int (B)\ (31) 

It depends on the flow parameter B and is anti-hermitian. The spit-up between 
the noninter acting and the interacting part of the Hamiltonian has to be 
defined throughout the flow. For the Hubbard Hamiltonian, this is simply 
achieved by promoting the parameters in the Hamiltonian to 'flowing', in- 
dependent variables. Hence we start with the following ansatz for the flowing 
Hubbard Hamiltonian 



H(t>0;B )= e k (B) :c\ a c ka :+ u wi2{ B ) rc^d-rcj,,^: S 1+2 



l'+2' 

u\i2,'V2\J-> ) •^i/-|-^lt L; 2 / | t ' 2 J-' u 

fceK,cr6{T,l} 1'12'2SK 



(32) 



Since the flow of the Hamiltonian is closely related to the net energy differ- 
ence of scattering processes, an iterative approach will show that the flowing 
interaction strength inevitably depends on Hamiltonian energies Ux>i2>2(B) = 
U(ei', e 2 >, ei, e 2 ; B) with an initial condition £/i>2'i2(A) = 0) = U. Inserting (32) 
into (31) makes the canonical generator more explicit. With Aei'i2'2(-S) = 
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e v (B) - ei(B) + e 2 >(B) - e 2 (B) it reads 



77(B) = Yl Uvw 2 {B) Ae V i 2 , 2 (B) :c\,^c{, i c 2l :5\+ 2 2 ' (33) 

l'12'2eK 

Still this is an implicit definition of the generator since the functional form 
of the flowing interaction strength is not known explicitly. Hence an iterative 
approach to the correct and consistent definition of the generator is necessary. 



It starts with a first parametrization of the flowing coupling constants in ( 33 ) . 
In a second step its explicit action onto the Hamiltonian allows to calculate 
an improved parametrizatior p^] 



5.1.3 Differential flow equations 

Let U(B, dB) = e r, ^ B ' dB be an infinitesimal unitary transformation promoting 
the Hamiltonian or any other quantum mechanical observable O(B) to a new 
representation 0(B+dB). A leading order expansion of its action in dB (which 



corresponds to the angle <p in (2.2.1)) 

0{B + dB) = U(B, dB)0(B)U-\B, dB) « 0{B) + [rj(B), 0(B)] dB 
gives rise to a differential flow equation 

dO{B) 



dB 



lv(B),0(B)) (34) 



As a differential statement this is always exact and corresponds to the generic 
definition of a 'generator' in the theory of unitary operations. Note that in 



(2.2.1) second order contributions in the angle ip have been considered which 
have no equivalence in a differential approach. Approximations enter via the 
iterative interplay between the generator, the Hamiltonian and other observ- 
ables. The flow equation method is exact as long as they are considered as 
abstract, implicitly defined objects. Represented in a truncated multi-particle 



operator basis the commutator on the right hand side of equation (34) may 
generate new terms which have not been part of the original representation. 
Although the later may be extended by these operator terms, this typically 
runs into an infinite regression and requires approximate truncations. If such 
a truncation scheme has been established the differential flow equation for 
operators can be decomposed into a set of coupled differential equations for 
flowing scalar parameters. We show this first for the transformation of the 



13 This interplay between the definition of the transformation and its action onto 
the Hamiltonian is a characteristic trait of the flow equation technique. It comes as 
a consequence both of the demand that the transformation should diagonalize the 



Hamiltonian and of the very generic construction of the canonical generator in (31 ). 
The implicit definition allows for a generator which is intrinsically well-adapted to 
the particular Hamiltonian. 
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Hamiltonian which serves to make the canonical generator explicit. All ap- 
proximations are made with respect to a perturbative expansion of the initial 
interaction strength U = U(B = 0). 



5.1.4 Transformation of the Hamiltonian and the generator in leading order 

We start with the Hamiltonian J32|) a nd a straightforwardly parametrized 
generator rj^(B) = 77(0) oc U [c/p3J)]. The first order contribution to the 
flow equation for the Hamiltonian, i.e. with 0(B) = H(B), comes from the 
commutator [rj(0),H } oc U and reads 

dUvw2(B) _ , , 2 

— — U (Aei/i2'2j 

It can be integrated and gives a leading order parametrization for the depen- 
dence of the flowing interaction strength on energy and on the flow parameter: 

U VV2f2 {B) = U e -(A^ 12 , 2 ) 2 B (35) 

With this parametrization the first-order approximation of the canonical gen- 
erator reads 

V W(B) = U £ Ae vlV2 e~^'^ 2B x 4^4^: (36) 

l'12'2eK 

Its characteristic feature is the exponential cutoff function. It suppresses in- 
elastic scattering processes which violate energy conservation on an energy 
scale set by ~ B~ l l 2 . 

Higher order contributions to the flow of the Hamiltonian include, for example, 
the renormalization of the one-particle energies or second order contributions 
to the flow of the interaction. We will show later that within the accuracy of 
a second order calculation of the time dependent number operator a leading 
order implementation of the diagonalization is sufficient; hence they can be 
neglected and we write (without any further calculation) the energy- diagonal 
Hamiltonian approximately as 

H(t>0;B = 00) = e fc :c Lcfe CT : + 

fceK )( re{T4} 

+ U ^( Ae i'i2'2)^ 1 1 ' + + fc fc 2 2 ' :c\^c 1] ci, l c 2l : (37) 

l'12'2eK 
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5.1.5 Flow equations for the creation operator 



In a second step we map the creation operator into the approximate energy 
eigenbasis of the Hamiltonian by means of the flow equation transformation 
(34 36). We define it as a basis independent observable CL which has different 



representations depending on the value of the flow parameter. It symbolizes the 
creation of a physical fermion. In the initial basis of noninteracting fermions, 
i.e. for B — 0, it coincides with the formal creation operator which we 
treat as the building block of an invariant many-particle operator basis. Rep- 
resented in this basis, CL acquires a composite multiparticle structure under 



the flow. New terms emerge on the right hand side of its flow equation (34) 
and mirror the dressing of an original electron by electron-hole excitations 
due to interaction effects. Respecting momentum and spin conservation, this 
motivates the ansatz 



Cj T (B) = h k ,(B) 4 T S 
+ E Mv2>i m {B) O :cI, T 4, iCu : (38a) 



l'2'l 



+ E M V2 , im (B) :c\,^c lV . (38b) 

l'2'l 

Here h(B) and M(B) are flowing parameters of the observable CL(S). The 

first can be linked to the quasiparticle residue via hk F (B) = ^Z(B) and de- 
picts the coherent overlap of the physical fermion with the related interaction- 
free momentum mode (quasiparticle) of the current representation; the later 
represents the incoherent background in the spectral function of an interacting 



system. Inserting this ansatz into (34) we calculate the differential flow equa- 



tions for the flowing parameters h(B) and M(B). A consistent normal ordering 
of all newly generated terms is essential and causes the emergence of character- 
istic fermionic phase space factors like Q122' = n\n 2 (l— n 2 ') + (l— %)(!— n 2 )n 2 i. 



= U E Ae fcl2 , 2 e -*<AW Q 122 , M 122 , m (B) (39) 
OB 12 , 2 

dM 5 ^(B) = _ u ^ _ n{2)] e -B { Ae^, &? M&22im{B) (4Q) 
Or) 2 , 2 

dMw^B) ^ h ^ B) A£5/56/i e _ B(Ae5 , 56 , l)2 (41) 



+ U y £ Hi)-n(2')} M 1&vm (B) Ae 2 , W5 e- B ^'^ 2 

12' 

12 

+ ^E [ M i5i'TTT - M 5'ii'TTt] K 1 ') - "(1)] Ae V i6'B e^**™^ 



i'l 
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A perturbative expansion in U of the flowing parameters h(B) and M(B) 
allows to reduce the complexity of the differential equations but depends on 
their initial conditions. Since the differential equations are linear, a solution for 
general initial conditions can be achieved as a linear superposition of solutions 
for independent initial configurations. We will discuss two cases: 

(A) Fully coherent initialization of one fermion in the momentum mode k: 
hii(B ) = Sf and M(B ) = for all possible indices 

(B) Fully incoherent initialization in the dressing state p'q'p Tll : hJB ) = 0, 
Mi'2'iih(Bq) = dl/dySi and M^(B ) = for all possible indices 14 



5.1.6 Case A: Perturbative analysis at the onset of the flow 



We discuss iteratively the action of the differential flow equations at the onset 
of the flow. In a first step, the flowing parameters h and M on the right 
hand side of the differential equations can be parametrized by their initial 



conditions (A). Hence only the first term at the right hand side of (41 ) remains 
influential and is, due to its pre-factor, of order U. Consequently, M<m(B > 0) 



is generated in first order of U. Re- inserted into (39) it accounts for a second 



order correction to the flowing parameter h(B). This describes the leading 



changes to the quasiparticle residue. Re-inserted into (40) and (41) it unfolds 



second order effects on flowing parameters of the M type. Since those do not 
influence second order results on the momentum distribution function (which 
will be shown later), this calculation can be, fortunately, dropped. Moreover, 
the initial conditions (A) are the natural ones to study the behavior of physical 
fermions. As they do not generate in relevant order, the ansatz for CL 



can be restricted to (38a). This is a consequence of the Pauli principle which 



disadvantages dressing of a fermion with excitations in the same spin state. 
Simplified flow equations read 



dh k ,(B) 
dB 



U 



E 

12'2 



Ae 



fei2'2 e 



-B(Ae fel2 , 2 )2 



Q122' M : 



122'TU 



(42a) 



OM ™*M a > = 17 £ * lt Ae, 5&1 c -sCAW (42b) 

We stress that this analysis is based on an approximate parametrization of 
the flowing parameters. It requires that their magnitude under the flow is 
still sufficiently described by the magnitude of their initial conditions. The 
following study will show that for the initial conditions of a physical electron 

14 To avoid confusion we stress that in the full context of our problem, this initial 
condition comes with a pre-factor proportional to U. Perturbative arguments always 
include it. 
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at the Fermi energy this remains the case throughout the flow and that, then, 
non-perturbative effects do not arise. 



5.1.7 Approximate analytical solution of the flow equations 



We integrate the simplified flow equations (42) with respect to the flow pa- 
rameter B. Similar to the above analysis we iteratively develop the flowing 
parameters from their initial conditions (A). We apply the same parametriza- 
tion h\-\ = 5i and integrate (42b) from Bq = to Bf = B. 



1 _ e -B(Ae 5 , 56 , k f 



Ae 5 , 



56'fc 



(43) 



This serves as an improved parametrization in (42a) and allows to write down 
a formal second order correction of h^. 



h (A) 



U 2 



(1 - p -B(e k -E) 2 \2 



(44) 



with a phase space factor 



I k (E) := ]T Q V2 >i6(E - e v - e 2 > + e^l 



+2' 
+k 



(45) 



l'2'l 



5.1.8 Evaluation of the phase space factor 



For a further understanding of the flow of our focus is on the energy 



dependence of the phase space factor Ik(E). Momentum conservation allows 
to eliminate the momentum index 1 = I' + 2' — k; we keep it as a shorthand 
notation. We evaluate Ik(E) in energy space assuming a constant density of 
states p 



/oo roo 
de v \ de 2 >Q(e r , e 2 >, e x )8{E - e v - e 2 > + ei) 
-oo J —oo 



(46) 



and observe that, at zero temperature, the phase space factor Q(ei,e2',e2) 
vanishes for all but two configurations of its energy arguments [69J: 



Case 


ev 


e 2 < 


ei 


Q 


E — By — 6\ + 6 2 i 


(a) 


> £f 


> e F 


< e F 


l 


E = + ei'| + |e 2 '| + |ei > 


(b) 


< £f 


< e F 


> e F 


l 


E = — ei/ — \e 2 t\ — ei < 



In both cases the limits of the energy integrations in (46) can be restricted by 



an approximate evaluation of the delta function. Its argument can be rewritten 
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by considering the signs of the one-particle energies e« which is done in the 
last column of the above table. From there one reads off that, since E is a 
sum of either solely positive or solely negative summands, E forms an upper 
or lower bound on both energy integrations individually. Then both cases (a) 
and (b) lead to the same approximate expression for the phase space factor 

I k (E) = p 2 [ E ' de v de 2 , 1 = p\E - e F f (47) 

We finally remark that the phase space factor Q{ey, e^-, ei) mirrors the particle- 
hole symmetry of the Hubbard model. It implies that all odd powers of a 



perturbative expansion vanish [64J. We will see in (5.3.5) that Q{ey, £2', £i) is 



in particular, responsible for the suppression of secular terms. 



5.1.9 Analysis of the later flow of h k ^(B) 



The approximation (47) allows to interpret the physical relevance of the for- 
mal result (44): At the Fermi surface, i.e. for k = fop, and at zero temperature 



its quadratic divergence is cancelled exactly by the phase space factor. This 
implies that the corrections to /ij^l remain small in second order of U through- 
out the flow. Hence the parametrization h k A \ = 1 used in (5.1.6) and (5.1.7) 



is justified even beyond the onset of the flow and the approximate analytical 
solution gives a correct description of the changes to the quasiparticle residue 
for all values of B. 



Away from the Fermi surface, however, there is no such cancellation by a phase 
space factor. Now the energy denominator indicates nonperturbative effects. 



Expanding the exponential in (44) shows that the numerator smoothes the 



peaked energy structure on a -B-dependent scale such that a pole only emerges 
in the limit of infinite B. Nonetheless, the contribution of an environment 
around e^kp to the energy integral in (44) grows under the flow until it infers 



a nonperturbative correction to h^ kF ^. This implies a full decay of a particle 
into incoherent excitations. The corresponding scales can be extracted from a 
first order expansion in B\ 15 | Since this is only justified for B{e k — E) 2 < 1 
we continuously restrict the environment [e^ — 1/ \/B, e& + 1/ \fB\ around the 
later pole. With e k = E — e& and (47), the correction term in (44) has the 



5 A similar analysis based on scaling out the independence from the energy integral 
leads to the same scales but requires a discussion of logarithmic and power-law 
divergences. 
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following dominant behavior for large values of B 



Ah (A) 



t(*) 



-XJ 2 B 2p - 



de k (e k + e k - e F f e\ 



2^3 



-U z p 



5v^B 3 1 



(48) 



The first term can be neglected for large B while the the second term exhibits 
the decay scale 



B* 



U*p«(e h - e F ) 



or E* = U z p'\e k - e F 



(49) 



The quadratic energy dependence of E* resembles the characteristic finite life- 
time of Landau quasiparticles away from the Fermi surface r ~ (e k — t F )~ 2 . 
Here it illustrates the complete decay of a physical fermion. This implies that 
on a scale set by B* the parametrization of = 5^ in (5.1.6) and (5.1.7) 



becomes unjustified and the analytical solution (44) breaks down. Instead 



a better solution of the differential flow equations would show a complete 
transfer of the spectral weight from the one-particle to many-particle repre- 
sentations. 



This allows to discuss the reliability of our solution: For particles in a pertur- 
batively small environment around the Fermi energy, i.e. e k — t kp ~ U, the 
energy scale E* ~ 0(U A ) is beyond the resolution of a second order calcu- 
lation. Hence the particle decay cannot be observed on all accessible times. 



Therefore the extension of the approximate solution (5.1.7) to an environment 



around the Fermi energy of radius U leads to a consistent second order result. 
Equilibrium calculations have shown that the broadening of the momentum 
distribution due to interaction effects is concentrated in this region [6Tf7T] . 



Outside of this environment, particles decay completely into incoherent mul- 



tiparticle excitations. This implies that, strictly speaking, the validity of (43) 
breaks down. Nonetheless, it shows that under the flow the spectral width 
Fm(B) of newly arising incoherent M-terms continuously decreases. This in- 
dicates that at the decay scale spectral weight is transferred from coherent 
particles (with a sharp spectral distribution) primarily to excitations with 
Tm ~ E*. Since even far away from the Fermi energy this is perturbatively 
small in U 2 , the effects of this small widening will not influence the shape of 



the momentum distribution function. Therefore the solution (5.1.7) allows for 



a calculation of the momentum distribution function on all energies. 
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5.1.10 Composition of the number operator 



Since the transformation for the creation operator of an electron around the 
Fermi energy has been established in second order perturbation theor y, th e 
number operator A4 = C k C k can be easily composed from the ansatz (38a). 
The momentum distribution is given by its expectation value with respect to 
the equilibrium ground state or the nonequilibrium initial state. Then only 
quadratic terms in h(B) and M(B) contribute 1 ' 



N k (B,t) = (FG\Af k (B,t)\FG) = 

= \h k (B,t)\ 2 n k + J2 n v n 2/ (l- ni ) \M* 2 , im (B, t)f S]^+ k 2 ' (50) 

l'2'l 

A possible time dependence has already been included for later reference. 



5.2 Equilibrium momentum distribution 



Similarly to the approach in [72] we observe that the flow equation transforma- 
tion resembles in many aspects a unitary implementation of Landau's theory of 
a Fermi liquid. Although a strict identification of the diagonal degrees of free- 
dom obtained from the flow equation approach with Landau quasiparticles 
is not possible, they motivate an analogous picture of flow equation quasi- 
particles; like their Landau counterparts they are stable only at the Fermi 
energy and subject to a residual quasiparticle interaction which is carried by 
the nonvanishing two-particle component of the energy-diagonal Hamiltonian. 
Nonetheless, these quasiparticles absorb most of the interaction effects into 
their definition such that the quasiparticle representation of the interacting 
ground state is, in good approximation, the filled Fermi sea. 



Hence the momentum distribution in equilibrium can be calculated in a similar 
way as the equilibrium occupation of the squeezed oscillator (13) using (50), 
(|43]) and (ph. 



N^ QU = (Q\ Cl(B = 0)C k {B = 0) \Q) = (Q \C k {B -> oo)^C k {B -> oo) |O ) = 

--Wj-iEWrt (51) 
J-oo [B - e k ) z 



16 This simple representation of N k (B, t) depends in a subtle way on the chosen nor- 
mal ordering prescription. Here the reference state of the normal ordering procedure 
is identical with that state for which the expectation value of the number operator 
is taken. If both states disagree, cross-terms proportional to ~ h X M appear. 



35 



A 



B 




A 




C 




K 




W 




A 




R 




□ 





B 


M» 


A 




C 




K 




W 




A 




R 




D 





Fig. 3. In nonequilibrium, a fermionic particle C\. decays with time. We only con- 
sider processes which contribute to the momentum distribution up to second order 
in U. In this schematic diagram the sequence of transformations applied to ct 
is illustrated. The FORWARD transformation starts with the fully coherent initial 
conditions of a physical particle. TIME arrows indicate the insertion of phase shifts, 
the BACKWARD transformation can be decomposed into fully coherent (A) initial 
conditions and fully incoherent (B) ones. The result of the sequence of transforma- 
tions is expressed in total coefficients h and M. 

where the phase space factor 

J k (E; n)=J2 & 5(e r + e v -e 1 -E)x 

l'2'l 

x [n k nx(l - ni')(l - n 2 >) - (1 - n k )(l - n 2 )n v n 2 >] (52) 
resembles the quasiparticle collision integral of a quantum Boltzmann equation 



. We will compare (51) with the time dependent momentum distribution 



after an interaction quench. 



5.3 Nonequilibrium momentum distribution 



We continue our approach towards the nonequilibrium momentum distribution 
by evaluating the time dependence of the creation operator represented in the 
initial basis (B = 0) 



0,t) = M0,t)4 T + E M i'2'im(0,t)O :c 



l'2'l 



t J 



according to fig. [I] Our aim is to write these time dependent coefficients h(Q, t) 
and M(0, t) in section (5.3.3) as functions of analogously defined (but time 



independent) parameters of the forward transformation (FT; see |5.1.7) and 



5.3.2 


) as 


ion ( 


5.3.1 



the used notation and the power counting of the perturbative analysis. 
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5.3.1 Time evolution 



We time evolve the creation operator in its energy- diagonal representation 
Ck]{B — > oo) with respect to the energy- diagonal Hamiltonian. Its leading 
part, the noninteracting Hamiltonian, generates a time evolution according to 
K(t,t Q = 0) = e iHot which can be treated exactly. Corrections arise due to the 
energy-diagonal interacting part of the Hamiltonian and cause the appearance 



of secular terms. It will be shown in section (5.3.5) that they can be neglected 



in a second order calculation. Hence the action of U onto the ansatz (38a) 



C^(B,t) = h[f FT {B)U\t^U{t^) 

+ M $fm( B ) « f (*, 0) :c\,^ 2ll c n :U(t, 0) (53) 



l'2'l 



is for B — > oo fully described by additional phases accompanying the param- 
eters hW' FT and M^ FT . 

h k1 (B -»• oo, t) = e itkt h$ ,FT (B -> oo) (54a) 
M V7fim {B -»• oo, t) = e^ e i' +£ 2'- £1 )* M&^fl -> oo) (54b) 

5.3.i? Inverse transformation 

The final step depicted in fig. [T] is the backward mapping of the time-evolved 
observables into the original representation of physical fermions. It is imple- 
mented by the inverse sequence of differential unitary transformations and 
simply achieved by interchanging the limits of the ^-integration in the for- 
ward transformation or, intuitively, by 'running the transformation backwards' 
|4Tj . Yet the decay of a physical fermion under the forward transformation has 
generated nonvanishing incoherent contributions in order U; therefore differ- 
ent initial conditions for the backward transformation of C^{B = oo,t > 0) 



apply. According to (5.1.5) a linear combination of the solutions for (A) and 



(B) leads to a full sol ution at arbitrary initial conditions. Case (A) can di- 



rectly be taken from (5.1.7) with M^' BT = —M^ AS>,FT and invariant h^ A \ 
The discussion of case (B) is simplified by noting that the pre-factor of its 
nonvanishing initial condition is proportional to U because of the generation 
of this term under the forward transformation. We will consider this additional 
power of U in a perturbative analysis of the relevant contributions but solve 
the differential equations for the unweighted initial conditions of case (B). 



Case B: Perturbative analysis and approximate solution 

At the onset of the backward flow, i.e. for large values of the flow parameter 
(B = oo), we insert the weighted initial conditions of case (B) as a constant 
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parametrization into the right hand side of the flow equations (39 41). With 
Mfjj ~ 0{U), h k ^(B) is generated by (39) to second order in U ; hence, within 
a second order calculation, the parametrization = holds throughout the 
backward transformation. Consequently, corrections to M^n are second order 
in U, as well is the generation of M^. Looking at (56b) and back to (50) 



shows that a second order result of the momentum distribution only requires 
the knowledge of M to first order. This is a priori known by the weighted 



initial condition. Integrating (39) gives 



h (B),BT 
,L k} 



M (B),BT (B 



(B) — U Qp'g'p 
V°2'°l 



e - B ( Ae p'pq' 

Atp'pq'k 



(55a) 
(55b) 



5.3.3 Composite transformation 

We finish the computation of the time dependent creation operator by com- 
posing the forward transformation (FT), the approximate time evolution and 
the backward transformation (BT) and represent the joint result in terms of 
time dependent parameters h k ^(B = 0, t) and Mm(B = 0,t). Fig. [3] gives a 
pictorial representation of these expressions. 

M* = 0,t) = h^e^h^+Y,^e^^M^ (56a) 

p'q'p 



M?, vim (B = 0,t) 



M (A),BT iekt 7 (A), FT V M {B) ' BT 
iW l'2'lTU e n k\ + IV1 p'q'pUl 

p'q'p 



r+ej-t 



m v'q'pni 
(56b) 



5.3.4 Nonequilibrium momentum distribution 



Inserting the results (43, 44, 54 and 55) via (56) into (50) yields the nonequi- 
librium time-dependent momentum distribution function for the initial state 

\n ). 



N? EQ (t) := (n \AT k (B = 0,t) |O ) = n k -AU 2 



*> sin 2 ((E - e fc )t/2) 

(57) 



We define the correlation-induced time-dependent correction to the momen- 
tum distribution 



AN™ EQ (t) = N™ Q (t)-n k 
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and perform a time average of its long-time limit. Note that equation (57) 
resembles the structure of Fermi's golden rule. There the same energy kernel 
consisting of the same sinusoidal time dependence and energy denominator 
becomes increasingly localized with time such that, in the limit of infinite 
time, the energy integral collapses into a secular term oc t [Hj. Here, however, 
the energy dependence of the phase space factor compensates for the energy 
denominator at the Fermi energy. This is a particular feature of a many- 
body fermionic system at zero temperature. Interchanging time average and 
energy integration then results in a factor of 1/2. Around the Fermi energy, the 
long-time limit of the second order perturbative correction to the correlated 
momentum distribution is given by 

AiV»:= hm(AiV fc NEQ W)* = -2f/ 2 1°° f^'-l ^ 2 AN^ QU (58) 

*^°o " ' J-oo [€k — E) Z 

Compared with the equilibrium result one observes the key factor two. It rep- 
resents a similar mismatch of the quasiparticle residue: Its correlation induced 
reduction is doubled in nonequilibrium compared to the equilibrium result 
1 _ z NE Q = 2(1 - Z E ^). 



5. 3. 5 Vanishing influence of leading secular terms 

The interpretation of a perturbatively calculated long-time limit result must 
be based on an analysis of possible secular corrections. Such terms may arise 
from a simultaneous expansion in both the interaction and time as it was done 
naively by the zeroth order approximation Wb^oo(^,^o = 0) = e iH ( B -* 00 ) t 



e lH ° l in (5.3.1); they may render a calculation unreliable on short time scales. 



We, firstly, notice that up to second order the decomposition 



& iHt _ e iHot e iH int t e t 2 [Ho,H int ]/2 ^ e iH t e iH int t 



is exact since in the energy-diagonal representation the generator rj(B) 



[H ,H int ] [cf. (33)] vanishes. 



Secondly, we will show that the first order expansion of AU(t,t = 0) 



e iH int (B->oo)t ^ ac ti n g on to the ansatz (38a) like Li in (53), does not influence the 



momentum distribution. The further analysis is simplified by the observation 
that these correction terms (which, of course, solely root in the time evolution 
of the creation operator) can be formally written as time-dependent correc- 
tions to the forward transformation. This allows to straightforwardly evaluate 
their contribution to the full time evolution of the momentum distribution 
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function by inserting the following corrections into (56) and (50): 



(13 



-it U M pVp FT( 5p'P9' 5 p'+V( Ae P'P9'A 



p' q'p 



it u z 



dE^^6(E 
E - e k 

AM^j = -it U h$6(Ae 1 wk)6l? k 



+2' 



(59a) 
(59b) 



(59a) vanishes because of the phase space evaluation presented in (5.1.8). The 
same argument holds for the influence of (59b) either directly in (56a) or 
finally in (50). Altogether, there is no leading correction. 

Moreover, there is no second order secular term contributing to the momentum 
distribution. The only one conceivable can be an outcome of the following 
transformation 



c{(B = 0) ^ c{(B = oo) 44, Cl : i D x ) 



HixA 

0(Ut) 



oo) ?5 c\{B = 0) 
o(i) 



The corresponding correction can be constructed from (59a) and (59b) and 
vanishes similar to (59a). 

Here we observe the suppression of secular terms due to the interplay of time 
evolution with respect to an energy diagonal Hamiltonian and fermionic phase 
space factors. It occurs, in most cases, already on the level of the transforma- 
tion of the creation operator. This illustrates the advantages of the chosen 
transformation scheme. One secular term, however, only vanishes because of 
the particular structure of the transformed number operator. 

We conclude that the time evaluation with respect to the noninteracting 
Hamiltonian Hq in (5.3.1) is justified and that the second order long-time 
limit (58) is not modified by secular corrections. 



5.4 Nonequilibrium energy relaxation 



Before we discuss the time behavior of the quenched Fermi liquid we establish 
some results on momentum mode averaged energies like the kinetic energy or 
the interaction energy. 

E^^(t):=U(Q \H int (t)\n ) (60) 

/oo 
de k e k N k (t) (61) 
-00 

Although the total energy of a closed system is conserved, its partition onto 
interaction energy and kinetic energy after the quench is time dependent. 
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Fig. 4. Sketch of the total energy (full line), kinetic energy (broken line) and in- 
teraction energy (dotted line) for the equilibrium (thin lines) and nonequilibrium 
(thick lines) case in second order perturbation theory. Nonequilibrium energies are 
compared at the quenching time and in the limit of infinite time. The bowed arrows 
indicate the corresponding energy relaxation, the straight arrows the gain in kinetic 
energy (broken filling) which equals the total excitation energy over the equilibrium 
ground state (full filling). 

The energy zero point is denned by the Fermi energy of the noninteracting 
Fermi gas ep = (Qq \ H \Q ) = 0. It agrees with the kinetic energy shortly 
after the quench E NECi ' Km (t = 0+) = because the state of the system 
is not changed by the quench directly and the time evolution of H with 
respect to H has not been effective. Moreover, as H mt is normal ordered with 
respect to \Qq), also the interaction energy shortly after the quench vanishes 
£NEQ,iNT( t = 0+ ) = o. This implies that the total energy before and after the 
quench remains identically zero _E NE Q = (Q Q \ H \Q ) := 0. However, since the 
ground state of the interacting normal-ordered Hamiltonian is lower that that 
of the Fermi gas, the system is initialized in an excited state. 



5.4-1 Equilibrium energies 

The point of reference of all energy considerations at nonnegative times, how- 
ever, is the total energy of the equilibrium ground state of the interacting 
system. We calculate it in perturbation theory with respect to the noninter- 
acting Hamiltonian H ; note that its eigenstates \m) (m G Nq with \Q ) = |0)) 
are non-degenerate. 

E E ^ := (<l\ H \Q) P = T - (0| tf |0> + (0| H int |0> + £ \Ml^Ml + (U 3 ) 

' -^T o ' rnJO ^ ^ f 



_£EXC 



(62) 



This allows to read off the excitation energy of the quenched quantum system 
above the interacting ground state E EXC := E NEQ - E EQV = U 2 p F a > 0. 
It is positive, second order in U and its precise value depends - due to the 
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summation over all quantum numbers m - on the lattice structure. We hide 
such details in a numerical prefactor a. 



Further equilibrium energies can be calculated from the Feynman-Hellman 
theorem [75 ] 17 1 

8E *T (U) = <n ( g)| gg^ mm) B 

rr pEQU.INT 

= (n{u)\-^P{u)) = —fj— (63) 

With (g and (g we know ^equ.int = _ 2 £Exc and £ E qu,kin = ^exc 



5.4-2 Nonequilibrium energies 

Next we calculate the nonequilibrium energies in the long-time limit of the 
second order calculation. As the Fermi energy is set to zero, 

£NEQ,KIN = j dtk tk XaTNEQ m 2jB EQU,KIN = ^EXC (g4) 



and £JNEq,int _ —2E EXC . All energies are sketched in fig. |4j Hence a second 
order calculation shows that the excitation energy of the quench is fully con- 
verted into add itional kinetic energy such that E'Neq.kin = ^equ.kin + ^exc 
while _E ne Q,int = ^;equ,int_ ^his j s a remarkable observation. 



5-4-3 Interpretation 

The redistribution of energy between potential and kinetic energy occurs in an 
early phase of the evolution of the model. Although the full excitation energy 
is transferred to the kinetic degrees of freedom this does not mean that the 
model has thermalized. While on average the kinetic energy has acquired its 
final value the distribution of the energy on the various (momentum) degrees 
of freedom is still a nonequilibrium one. This can be seen in the momentum 
distribution function which will relax to thermal equilibrium only on a much 
longer time scale. In between there is a time regime in which a characteristic 
nonequilibrium momentum distribution is retained while the energies have 
already relaxed. 



17 



Note that the norm of the interacting ground state is invariant. 
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Fig. 5. The time evolution of the momentum distribution A NE Q(e) is plotted around 
the Fermi energy for ppU = 0.6. A fast reduction of the discontinuity and 1/t-oscil- 
lations can be observed. The arrow in the plot for t = 2.5 indicates the size of the 
quasiparticle residue in the quasi-steady regime. 

5.5 Discussion of a multi-step dynamics of the nonequilibrium Hubbard model 




In the following we will illustrate the consequences of the above calculations 
and develop a picture of the real-time evolution of a Fermi liquid after an 
interaction quench. We observe three different time regimes which we illustrate 



by a numerical evaluation of the time dependent momentum distribution (57). 



For computational convenience we use the limit of infinite dimensions where 
momentum sums can be evaluated by energy integrals. It is generally assumed 
that in this limit the generic features of a Fermi liquid are retained. A Gaussian 
density of states p(e) = exp (— (e/t*) 2 /2) / \f2ixt* includes the constraints of a 
hyper cubic lattice, t* is linked to the hopping matrix element by dimensional 



scaling t — ► -4== to retain a nontrivial relation between the kinetic and the 

*/2ct 



interaction energy in all dimensions [76]. pp = p(e = 0) denotes the density of 
states at the Fermi level. For three time steps explicit results are depicted in 



fig. 5.5 



5.5.1 Short-time quasiparticle buildup and nonequilibrium state 

The first time regime is the one covered by the above second order perturba- 
tive calculation. The evolution of the time dependent momentum distribution 
for physical fermions exhibits the build-up of multi-particle correlations from 
a noninteracting Fermi gas which leads to the formation of a quasiparticle 
description. The shrinking of the Fermi surface discontinuity of the momen- 
tum distribution to a finite, nonzero value (Fig. 5.5 mirrors the fast reduction 
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Fig. 6. The universal curves for the correction to the momentum distribution 
ANk = Nk — rik are given for both the equilibrium and for the nonequilibrium 
quasi-steady state in the weak-coupling limit. 

of the quasiparticle residue from one to Z NE ® < 1 on a time scale is set by 
second order perturbation theory < t < p^ l U~ 2 . 1/t-oscillations accompany 
this process. 

Afterwards, the momentum distribution functions indicates that a zero tem- 
perature Fermi liquid description holds. Hence, quasiparticles are well-defined 
quantities around the Fermi surface. However, the mismatch of the quasipar- 



ticle residue by a factor of two (58) when compared to the corresponding 
equilibrium distribution displays its continued nonequilibrium nature and a 
postponed relaxation. 

It is helpful to change the point of view from a momentum distribution of 
physical fermions (PF) to one of quasiparticles (QP). They agree with each 
other on the main features. However, since interaction effects are absorbed into 
the definition of the quasiparticles, the quasiparticle distribution exhibits less 
pronounced correlation-induced signatures. In equilibrium they vanish com- 
pletely and a zero temperature quasiparticle distribution always equals a filled 
Fermi sea with a Fermi surface discontinuity of size one = 0(6^ — 6^). 

In nonequilibrium, this absorption of interaction effects is incomplete. Due to 
the increased correlation-dependent reduction of the quasiparticle residue by 
a factor of two a mapping of the nonequilibrium momentum distribution for 
physical fermions into a quasiparticle representation will generate a distribu- 
tion with a reduced Fermi surface discontinuity. In second order perturbation 
theory A^ E = iV^ u ' PF . The reduced discontinuity of the quasiparticle 
distribution now describes the deviation from equilibrium. We will discuss its 
further relaxation. 

While, so far, we have only observed a partial relaxation of the momentum 
distribution, the second order calculation shows a complete transfer of the ex- 
citation energy from interaction energy to kinetic energy [cf. fig. [2]. However, 
since the momentum distribution still exhibits zero-temperature features, this 
is not related to heating. The simultaneous relaxation of average energies and 
non-relaxation of other, mode-specific expectation values is known as prether- 
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malization. It has been described in nonequilibrium quantum field theories 
modeling, for instance, the early universe [77J and underlines the nonequilib- 
rium nature of the final state depicted by fig. |5.5 



5.5.2 Intermediate quasi-steady regime. 

So far we have observed the build-up of a characteristic nonequilibrium state 
which we have characterized by its energies and momentum distribution. A 
study of corrections to the second order perturbative result shows that it ex- 
tends throughout a second time regime: As the Hubbard model is particle-hole 
symmetric [64J the next nonvanishing contributions are expected in fourth 
order. Their relevance for the further dynamics will be discussed later. Here 
we note that there are no immediate changes to the state of the system for 
times t > p^U~ 2 . Hence the nonequilibrium Fermi liquid state represents 
an intermediate quasi-steady regime of the dynamics; for small values of the 
interaction U it holds for p^U^ 2 < t < pp 3 U~ 4 . This should simplify the 
observation of this nonequilibrium state in prospective experiments. 



5.5.3 Long-time behavior - Thermalization 

The mentioned fourth order corrections to the flow equations calculation de- 
scribe the characteristic dynamics of a third time regime on a scale t ~ Pp 3 U^ 4 . 
It originates both from higher order terms of the diagonalizing transformation 
and from the full time evolution generated by the energy diagonal interaction 
Hamiltonian. Their full calculation is beyond the scope of this work. Fortu- 
nately, an effective treatment of elastic scattering processes in a quasiparticle 
representation is provided by the quantum Boltzmann equation (QBE) [75] 

^p!) = _ pp u* Jk{E = 6k7 ivQP (t)) . (65 ) 

Based on its main prerequisite, a well-established quasiparticle picture, it 
models phenomenologically the relaxation of a nonequilibrium quasiparticle 
distribution to a thermal one. One generally expects that it describes the cor- 
rect long-time behavior although its proper derivation remains an unsolved 
but well-discussed problem of mathematical physics; this contains the delicate 
question how the transition from a deterministic quantum dynamics to a ir- 
reversible statistical description can be rigorously achieved [78]. Here we only 
motivate its application. 

The characteristic features of the quantum Boltzmann equation can be read 
off its right hand side which is commonly referred to as the scattering inte- 
gral |79J. Since Jk{E = e&,n) is energy conserving it describes similar elastic 
two-particle scattering processes as the energy-diagonal interacting part of the 
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Hamiltonian. Thus a (QBE) description appears as a natural extension of a 
flow equation analysis which only addresses energy non-diagonal processes. 
This motivation is backed by the observation of a related coincidence: a per- 
turbative expansion of the scattering integral for the intermediate steady state 
leads to an analogous fourth order correction as it is expected in the flow equa- 
tion approach. Therefore we link the QBE to the previous dynamics by taking 
the quasiparticle momentum distribution of the intermediate nonequilibrium 
Fermi liquid state jv^ eq ' qp as its initial condition. Because jyQ p - NE Q allows 
nonzero phase space for scattering processes in the vicinity of the Fermi sur- 
face, linearizing the phase space factor in the scattering integral shows that 
the initial quasiparticle distribution function starts to evolve on the time scale 
t oc pp 3 U~ 4 . This implies that the quasi-steady fermionic distribution function 



depicted in the last panel of fig. 5.5 starts to decay on the same time scale. 



The further dynamics of the quasiparticle momentum distribution function 
follows, again, from the scattering integral. Since Jk{E = €k,n) vanishes for 



Fermi-Dirac distributions (n = n FD ) these are the stable fixed points of (65). 
Moreover, the scattering integral conserves the kinetic energy such that the 
evolution towards a fixed point is constrained to an energy hypersurface in 
phase space. Hence, if the quantum Boltzmann dynamics continues until it 
reaches its stable fixed point, thermalization of the momentum distribution 
can be expected. This implies that the excitation energy, which has relaxed 
into an excess of kinetic energy already at an earlier stage, is redistributed 
among the momentum modes until a thermal distribution is achieved. The 
corresponding temperature T t ^ ~ U of the thermal momentum distribution 
follows directly from fitting its Sommerfeld expansion [SO] to the excitation 
energy. 

Notice that fourth order corrections to the diagonalizing transformation may 
cause an obliteration of the Fermi surface discontinuity even for short times. 
However, after the quench, the momentum distribution will still show a steep 
descent; therefore its widening can be safely neglected. We want to mention 
that the assumption of a quasiparticle picture for all later times is nontrivial. 
This is a general question of a QBE approach and is usually accepted. 



5.6 Consequences of the Hubbard dynamics 



In the past section we have observed the separation of two time scales of the 
Hubbard dynamics. While interaction effects lead to a rapid establishment of 
a quasiparticle picture, the equilibration of the momentum distribution, i.e. 
heating, is deferred to a much later time. This delayed relaxation is a con- 
sequence of two fundamental properties: Firstly, the Pauli principle imposes 
characteristic phase space restrictions on a multiparticle fermionic system. 
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Those suppress two-particle scattering processes and, thereby, reduce the effi- 
ciency of an interaction-driven momentum relaxation; they are responsible for 
the generic appearance of particle-like low energy excitation physics in a Fermi 
liquid such that a meaningful analogy between the behavior of the squeezed 
one-particle oscillator and the many-body Hubbard model can be justified. 
Secondly, translational invariance implies the conservation of lattice momenta. 
Hence there is no other way of momentum relaxation than by momentum 
transfer in two-particle (or higher) scattering processes. In combination, these 
two properties form a restrictive bottleneck for the relaxation dynamics. The 
particular form of the interaction, however, is less important. For the above 
calculations we have assumed the applicability of perturbation theory in the 
interaction strength. Since only second order and fourth order terms describe 
the evolution of the momentum distribution, there is no difference between 
attractive and repulsive interactions; moreover, a generalization to nonlocal 
interactions is easily possible by introducing momentum dependent interaction 
matrix elements. The main observation of a characteristic mismatch between 
the -interaction dependent- zero temperature correlated equilibrium ground 
state of the momentum distribution and a similarly shaped distribution in 
the intermediate regime of the nonequilibrium case persists. Hence we expect 
similar nonequilibrium behavior for a large class of weakly interacting and 
perturbatively approachable model systems independent of the exact nature 
and the particular form of the interaction. This reflects the rather generic ap- 
plicability of Fermi liquid theory for not too strongly interacting systems in 
equilibrium. 

Moreover, our findings are relevant for studies focussing on the nonequilib- 
rium physics of models with a Fermi liquid instability (FLI). Let us consider 
a quench from a noninteracting Fermi gas into a phase which exhibits such an 
instability. In the subsequent dynamics on a buildup time £fli-b one expects 
both the buildup of the Fermi liquid instability and characteristic nonequilib- 
rium physics related to it. Since non-perturbative weak interaction instabilities 
are typically linked to exponentially small energy scales, £fli-b will be large 
such that characteristic features of the instability are not observable for short 
times after the quench. In this regime, our perturbative calculation for the 
nonequilibrium Fermi liquid applies approximately even in the presence of a 
nonperturbative instability. Quenching into a FLI-phase then requires us to 
compare the timescales of the dynamics of the nonequilibrium Fermi liquid 
with that of the instability. 

If the FLI-phase is nonperturbative and distinguished by a gap in the energy 
spectrum, as it is, for example, the superconducting phase of a Hubbard model 
with an attractive interaction, an excitation beyond the energy gap is essential 
to observe any characteristic nonequilibrium behavior. This excitation can be 
induced by a sudden interaction quench. As we have seen the inserted energy 
causes heating effects which may wipe out all signatures of the FLI. Still the 
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delayed onset of heating in a Fermi liquid can open a time window for the 
observation of the nonequilibrium dynamics even in the FLI regime. 

A popular example for such behavior is the BCS instability. Recently the study 
of its nonequilibrium dynamics following a sudden quench in the BCS inter- 
action has attracted a lot of attention; depending on the precise conditions of 
the quench, for instance oscillatory behavior in the order parameter Abcs(^) 
has been found [81|82f22| . These studies only focus on the behavior of the 
(nonlocal) BCS Hamiltonian which is an effective low energy description of a 
superfluid. Since its dynamics is integrable a complete topological classifica- 
tion of the behavior of all excited states could be given [22J and no heating 
is observed. The actual experimental realization in optical lattices, however, 
only allows for a quench of the local two-particle Hubbard interaction. Aside 
from the emergence of an effective BCS interaction, the persistent influence 
of ordinary Fermi liquid behavior can be expected in such systems. Then a 
quench simultaneously initializes the nonequilibrium dynamics of the instabil- 
ity and heating effects. For a sudden quench heating dominates in agreement 
with [82J (since T e g 3> Abcs = exp(— 1/ and makes the nonequilibrium 

BCS dynamics unobservable. Hence a further analysis of the crossover be- 
tween instantaneous and adiabatic switching is required to study the visibility 
of such nonequilibrium dynamics. 



6 Conclusions 

In this paper we have presented the real-time dynamics following an interac- 
tion quench for systems with a discrete energy spectrum and for a Fermi liquid. 
In all cases we have discussed our key observation which is a discrepancy be- 
tween the equilibrium and the long-time averaged nonequilibrium occupation 
by a factor of two. This factor appears precisely in the modification of the 
noninteracting occupation due to interaction effects and illustrates a simple 
example for the interplay of interactions and nonequilibrium conditions. 

In a first analysis this factor was calculated for the squeezed oscillator. Since 
this is an exactly solvable one-particle model, a comparison between the per- 
turbative and the exact result showed that this discrepancy is not an artifact 
of perturbation theory. The precise value of two, however, is only reached in 
the limit of weak interaction. Afterwards, this observation was formulated as 
a theorem applicable to a larger class of discrete systems and observables. We 
have given two proofs which point out aspects of the origin of the factor two. 
One of them reconstructs and illustrates the applied transformation scheme 
which we used to calculate the nonequilibrium dynamics. 

These observations for discrete systems are paralleled by analogous behav- 
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ior for the zero-temperature momentum distribution function of a quenched 
Fermi liquid. In the second part of this paper we presented details on a weak 
interaction quench within the Fermi liquid phase of the Hubbard model in 
more than one dimension pp. To cope with the continuous energy spectrum 
of this Hamiltonian we based our calculation on the flow equation method 
following Wegner to achieve an approximate diagonalization in energy space. 
Since the Hubbard model in more than one dimension is nonintegrable, the 
expectation is that the expectation values of simple observables should ther- 
malize. Hence the characteristic mismatch of the nonequilibrium momentum 
distribution function when compared with the equilibrium one can only be a 
transient phase of a longer relaxation dynamics. Accordingly, we have found 
a three step dynamics of the quenched Fermi liquid: During the first phase a 
nonequilibrium quasiparticle description builds up which is characterized by 
an untypical value for the nonequilibrium quasiparticle residue. Again, the 
factor of two appears in the correlation induced reduction of the quasiparti- 
cle residue and mimics one-particle behavior. The momentum distribution of 
this transient nonequilibrium state remains frozen throughout a second phase 
which can be large for weak interaction. This has been explained by referring 
to the restricted phase space for two-particle scattering in a zero-temperature 
fermionic many-body system and corresponds to the collisionless regime. It 
delays momentum relaxation but does not prevent prethermalization of the 
kinetic energy. A third phase of the dynamics of the momentum distribution 
follows from the nonequilibrium nature of the transient state. Its nonequilib- 
rium momentum distribution opens phase space for further scattering events. 
Upholding a quasiparticle picture their backaction onto the momentum dis- 
tribution is described by a quantum Boltzmann equation which is an effective 
kinetic equation for the momentum distribution function in terms of a scat- 
tering integral. Since the later vanishes for equilibrium distribution functions 
and nonequilibrium derivations are small in U 2 it predicts the thermalization 
of the momentum distribution on a longer timescale set by £ t h ~ p 3 U~ 4 . 

This is a major difference to a similar analysis for fermions in one dimension 
[27] . There a quench in the forward scattering leads to an integrable dynamics 
which equals, after bosonization, that one of the squeezed oscillator for each 
of the momentum modes. In consequence, a similar mismatch between critical 
exponents of equilibrium and nonequilibrium correlation functions has been 
found which approaches the factor two in a perturbative limit. However, as no 
residual quasiparticle interaction occurs, thermalization is impossible: from 
our point of view our quasi-steady regime extends to all times in this one 
dimensional model. Comparing the results for the one dimensional case with 
our observations for higher dimensions elucidates how thermalization occurs 
or is inhibited in these two translationally invariant systems. 

Another perspective for further research lies in the examination of the crossover 
from instantaneous to adiabatic switching. 
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